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SUMMARY 


The general problem of heat transfer In a concentric 
annulus with hydrodynamlcally fully established laminar flow 
Is considered. Under the assumptions that the fluid in the 
annulus is Incompressible, has constant properties, and that 
internal viscous generation is negligible, the differential 
equation for the temperature field is solved for certain 
special or fundamental boundary conditions. These funda- 
mental solutions can be superposed to satisfy any general 
axially symmetric boundary conditions. 

The fundamental boundary conditions fall into four 
basic categories. The fluid has initially a uniform zero 
temperature; at some plane normal to the annular axis, 
x = 0, one of the following occurs: 

1. The temperature on one wall is increased to unity, 
the other wall is maintained at zero; 

2. The heat flux on one wall is increased to a constant, 
the other wall is insulated; 

3. The temperature on one wall is increased to unity, 
the other wall is insulated; 

4. The heat flux on one wall is increased to a constant, 
the other wall is maintained at zero temperature . 

Since each of these four cases is applied to both the 
inner and the outer walls of the annulus eight separate 
solutions to the energy equation are generated. From the 
fundamental solutions dimensionless wall temperatures, wall 
heat fluxes, and mean temperatures are defined and evaluated. 
The use of these quantities in solving problems of simple 
superposition is described and illustrated. 
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These fundamental solutions are obtained for values 
of the ratio of the Inner to the outer radius, r*, of 

0. 02, 0.05, 0.1, 0.25, and 0.5. These values together with 
existing results for r* of unity (the parallel plane channel) 
are presented in tabular and graphical form to facilitate 
interpolation so that problems in virtually any annular 
geometry can be solved. 

The solutions to the differential equation are obtained 
by the method of separation of variables. The wall boundary 
conditions are reduced to a homogeneous form by subtracting 
from the desired solutions the fully developed solutions, 

1. e., the solutions valid as the axial coordinate, x , tends 
to infinity. The eigenvalues, eigenfunctions, and expansion 
coefficients arising in the resulting Sturm-Liouville equation 
are obtained through a direct numerical integration facilitated 
by a convergent iteration method for obtaining the eigenvalues. 
The eigenvalues and relevant combinations of the expansion 
coefficients and the eigenfunctions are presented in tabular 
and graphical form for the previous values of r* and for 

r* = 0 (the circular tube). These quantities enable the 
consideration of continuous axial variation of boundary con- 
ditions via the Stieltjes integral. 

An apparatus for the experimental study of heat transfer 
to air with constant (but arbitrary level) heat flux at 
either or both of the boundary walls is described. The 
apparatus enables the direct experimental determination of 
the fundamental solutions satisfying boundary condition 2 
above, and superpositions thereof. 

A comparison of these experimentally determined values 
with those predicted by the theory is presented and 
exhibits the very good agreement between them. 
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I. INTRODUCTION 


I. A. The Problem 

Geometrically an annulus is a region bounded both in- 
ternally and externally by cylindrical surfaces. Complete 
geometric specification of the annulus requires not only 
the size of one of the tubes but the ratio of the radii 
and some knowledge of the orientation of the axes of the 
two tubes. In the concentric annulus the axes are coincident. 
Figure I.A.l illustrates a concentric annular passage and 
defines the coordinate system which will be used In descri- 
bing it. The radial qoordinate, r, Is measured from the 



common axis, the axial distance , x, is measured from an 
arbitrary reference plane normal to the tube axes. The 
concentric annular passage is a natural generalization 
of the circular tube and encompasses both the tube and the 
parallel plane channel as limiting geometries as the ratio 
of the inner radius to the outer radius, r*, varies between 
the limits zero and one. 

As a fundamental Internal flow geometry the annulus 
is frequently used In a situation where heat is transferred 
to or from the contained fluid. An annulus is a simple form 
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for a two fluid heat exchanger and a water or steam jacket 
around a pipe has been used for years as a method of heating 
or cooling the central fluid. An annular coolant passage sur- 
rounding a central fuel element has been used as a core con- 
figuration for a nuclear reactor. 

Since there are two distinct boundary surfaces annular 
heat transfer is inherently asymmetric, i.e., heat can be 
transferred at either the inner or the outer surfaces or both. 
In practice some heat transfer usually does occur at both 
boundaries even though it may be unwanted heat leakage through 
an insulating blanket or result from radiation from one surface 
to the other. Further, the temperature and/or the heat flux 
at either surface may vary both peripherally and axially in 
an arbitrary manner. 

The possible asymmetry in heating conditions and the 
relative geometric complexity combine to create a large num- 
ber of heat transfer problems in annular passages which are 
of significant interest. 

I. B. Method of Solution 

In order to formulate a systematic approach to the 
solution of problems in annular heat transfer we will restrict 
our considerations somewhat. Let us consider heat transfer 
In a concentric annular passage where the flow is hydrodynami- 
cally fully established and the temperature field is axially 
symmetric . 

The problems which might possibly be of Interest here 
can be defined by various conditions at the boundary surfaces . 
These are : 

1. Temperature arbitrarily specified on both surfaces, 

2. Heat flux arbitrarily specified on both surfaces, 

3- Heat flux arbitrarily specified on one surface and 
temperature arbitrarily specified on the other. 
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4. Heat flux specified over a portion of one or both 
surfaces and temperature specified over another 
portion. 

If we assume that the fluid is incompressible , has 
constant properties , and that axial conduction can be 
neglected, then the equation of the temperature field for 
laminar flow and no energy generated Internally is 

d 2 t t 1 dt _ U(r) dt ( I .B. l) 

dr 2 r dr a dx 

Since the energy equation is linear , any linear combination 
of solutions is again a solution. It may then, be possible 
to develop certain special or fundamental solutions to 
Eq. (I.B.1) satisfying particular boundary conditions which 
can be combined to satisfy the more general boundary condi- 
tions 1 through 4 above. 

For fully established, incompressible, turbulent flow, 
one can define an eddy conductivity, and the energy equation 
becomes 

|p [ r (c H + a) || ] = r U(r) || (l.B.2) 

Then if Is independent of the temperature, the equation 

Is still linear and all of the superposition techniques will 
be valid. Our consideration of the fundamental solutions is 
not restricted to laminar flow, in fact, Eq. (I.B.2) can be 
considered as including the case of laminar flow when = 0 

1 . The fundamental solutions of the first kind 

Consider the case where the fluid Is flowing with a 
constant temperature t g . At the station x = x q , the 
temperature on wall j (where J = i or o corresponding to 
the inner or the outer walls respectively )is increased to tj 
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and remains at that value for all x > x while the opposite 

wall is maintained at t . 

e 



FIG. I.B.l: The fundamental solutions of 

the first kind, temperature 
step at the inner wall. 

If we let 

= (I.B.3) 

*j - S 

where t is a solution to the energy equation satisfying 
the boundary conditions Indicated above, then 0^^ 
will be a solution of the energy equation satisfying the 
boundary conditions 

Sj 1 ) = 1 on wall j, 

= 0 on the opposite wall for x > x , 

(I.B.4) 


and 


8 (d . 

J 


= 0 for all 




at 


x = x 
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2. The fundamental solutions of the second kind 

Consider the case where the fluid again has the initial 
temperature t , and at x = the heat flux on wall j 

is increased to q^j , while the opposite wall is insulated. 


Insiilated 



PIG. I.B.2: The fundamental solutions 

of the second kind. Inner 
wall heated. 


Let 


e (2) 

J 


t " t 


q" D h /k 


(I.B.5) 


where = 2(r Q - r.^) , is the hydraulic diameter, and k is 

the fluid thermal conductivity. will be a solution of 

the energy equation satisfying the boundary conditions 

<o\ d 0 i 2) 

= _ n — — = 1 on wall j 

j h on 

= 0 on the opposite wall for x > 


and 


e< 2 ) 


0 for all r / Tj at x - x Q 


(I.B.6) 


Here n Is the direction normal to the wall, measured as 
positive into the fluid. 
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3 • The fundametal solutions of the third kind 


At x = x wall j is increased to 
o 

opposite wall is insulated. 


while the 



FIG. I.B.3: The fundamental solutions of 

the third kind., temperature step 
at the inner wall. 


Let 


S <3) = 

J 


t - t 


h - 


(1.B.7) 


9 


(3) 


will he a solution to the energy equation satisfying 


the boundary conditions 


*j 3) 

and 


= - D 


e<3) . 

h 3?T“ 


(3) 

0, 


0 for 


1 on wall j, 

= 0 on the opposite wall 

all r £ r . at x < x 
j — o 


for x > x , 
(I.B.8) 
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4 , The fundamental solutions of the fourth kind 


At x = x wall j has the applied heat flux q'j 
while the opposite wall is maintained at the upstream 
temperature t g . ^ 

fl* = 0 




FIG. I.B.4: The fundamental solutions of 

the fourth kind, inner wall 
heated. 


Let 


^ (4) 


t - t. 


*3 V 


k 


(I.B.9) 


0 3 (4) will he a solution of the energy equation 
satisfying the boundary conditions 

^ 4 > 




(4) _ 
J 


D 


h 2n 


J — = 1 on wall J, 


0^4) = o on the opposite wall for x > x Q , 
J 


and 6^ = 0 for all r £ r^ at x < x Q 


(i.B. 10) 


It is obvious from the defining equations that the fun 
damental solutions, 9^, are -» fact, dimensionless 
temperatures. Similarly the S>)j can be considered as 
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dimensionless heat fluxes. 

These four fundamental solutions can be superposed to 
form solutions satisfying the more general boundary conditions. 

The engineer is frequently not concerned with the details 
of the fluid temperature field but only with the temperatures 
and the heat fluxes at the boundary walls and their relation- 
ship to the mixed mean fluid temperature. It will be con- 
venient to introduce a generalized notation for the funda- 
mental solutions. The general non-dimensional wall temperature 
for any fundamental case, k, has the form efe) 

where 

& = 1 or o , referring to the wall at which the 
temperature is evaluated; 


j = i or o , referring to the wall at which the 
non-zero boundary condition is applied. 


Similarly the general non-dimensional wall flux has 
the form where the subscripts have the same meaning 

as before. The general mixed-mean (or cup mixing) temperature 
is denoted as where j = i or o referring to the 

wall with the non-zero boundary condition. By definition 


£>( f) 

mj 



/ U dA 

A J 


or 


9 


00 

mj 


(x) = 


tt Ir 


o 


- r?)U 


/ ° U(r) 0j k ^(r,x)2Trr dr 


m 


(I.B.11) 


where U m is the bulk mean fluid velocity. 
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We will, then, he principally concerned with constructing 
the functions giving the wall temperatures and heat fluxes, 
and the mean temperature in non-dimensional form for the 
various boundary conditions. We may summarize the principal 
functions of interest for each of the fundamental cases 
as follows: 

(l) First kind 


We seek 

9 h> 

mi 

*U) 

ii 

oi 


0 (D 

mo 

oo 

*U) 

io 


From the boundary conditions 




= et 1 ) 
00 

= 1 

and 


- 

oi 

= 0 

(2) Second kind 
We seek 

e<?> 

mi 

e (2) 

ii 

e (2) 

oi 


e (2) 

mo 

0 C2) 

OO 

a C2) 

io 

From the boundary 

conditions 



= *(2) 
OO 

= 1 

and 


- ®C2) 

io 

= 0 

(S') Third kind 
We seek 


e (3) 

oi 

$(3) 

ii 


e (3) 

mo 

a (3) 

io 

$(3 ) 

00 
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From the boundary conditions 



iV- 

e (3) 

00 

= i 


®oi ) = 

<ji(3 ) 
io 

= 0 

Fourth kind 



We seek 

0 (*) 
mi 


®(*0 

oi 


0 (“) 
mo 

eW 

00 

io 

From the 

boundary conditions 


*(*) = 
ii 

00 

= i 


II 

-=t* *H 
O 
CD 

«(4) 

io 

= 0 


To illustrate the use of the method of superposition 
and the utility of the pertinent non-dimensional functions 
let us consider as examples special cases of the arbitrary 
boundary conditions. 

1 • Constant but unequal wall temperatures 

We have the boundary conditions that for x > x 


t = t on the outer wall, 

t = t. on the inner wall, and 

t = "tg for all r at x < x q . It can be seen from the 
boundary conditions (l.B.4) that if we let 


t 


t e' 1 ) 

O O 


+ 


t e (1) 

Vi 


+ te 


(I.B.12) 


t will be a solution to the energy equation satisfying the 
desired boundary conditions. Since the wall temperatures 
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are prescribed we are principally interested in the heating 
rates at the walls and the axial variation of the mean tem- 
perature. By manipulation of (l.B.12) we can obtain 


and 


a" = - [ t + t ] 

q i D, i ii o 10 

h 


" = £ [ t. <&^ + t ] 


4 o D 


h 


i oi O' oo 


Because of the linearity of Eq. (i.B.ll) we have 


t = t J 1) . t 0 U) + t 

i mi o mo ( 


2 . Constant but unequal wall heat fluxes 
If we specify 


(I.B.13) 


(I.B.14) 


(I.B.15) 


q" = q'| on the inner wall, 

q" = q^ on the outer wall, for x > x Q and 
t = t for all r when x < x Q , then we can see from 
(I.B.6) that 


= IT [ q i 0 i 2) + Q " ] + 


L o o 


(I.B.16) 


is a solution of the energy equation satisfying the appropri- 
ate boundary conditions. Further we find 


D 

Sn = k 


h 


t q',' ^ + q„ 0 io } 1 + t 


D 

*i = k 


h 


l i mi 

[ q',' + q„ 1 + t 


1 11 


(I.B.17) 

(I.B.18) 


11 



(I.B.19) 


o 


D, 

= Jl [ 

k L 


q" 0 


( 2 ) 

oi 


+ q" 0 


( 2 ) 

oo 


] + t 


3* Consta nt heat flux on one wall and constant temperature 
on the other 

Let t = t on wall a, 

ci 

9. = Q-Vj on wall t> , for x > x and 

t = t e for a11 r when x ^ x Q . Combining the mixed cases 
we see that 


t 





+ t . 


(I.B.20) 


Is a solution satisfying the correct boundary conditions, 
the other quantities follow as before. 

Then If the j ^ ■> > anc ^ the wer> o available 

as functions of the axial coordinate, x , for the particular 
annulus of interest, the heat fluxes or temperatures at the 
walls and the mean temperature at any axial location could 
be computed readily by a simple linear combination of the 
appropriate values. 

We would like to generalize these simpler considerations 
to the case of arbitrary axial variation of the wall boundary 
conditions. Consider the wall temperature variation shown 
in Pig. I.B.5. One might visualize this shape as being con- 
structed from a large number of steps of infinitesimal height. 
Summing up the contributions of these steps and passing to the 
limit, the solution of the energy equation satisfying this 
boundary condition may be written in the form of a Stieltjes 
integral. 
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(I.B.21) 

The integral is evaluated by using dt^/d^ wherever t^(£) 
is continuous and adding 6 j ^ (r,x-£ n )At w ( £ n ) 
expression for t(r,x) for x > Bor a more detailed 

discussion of this superposition method see Klein and 

pp 16 

Tribus , or Reynolds et al. 

Using ( I.B.21) we can write for the cases of arbitrary 
wall temperature and heat flux the following equations: 


n 


PIG. I.B.p 


£=x dt ( ^ 

t(r,x) = / 0 1} (r,x -0 (• w 

1=0 J 


"dT 


)de 


£=x / >, dt . ( £) 

t(r,x) = f ey* (r,x -0 ( —)d£, 

i=0 J ds 


+ ^ # (D (l , x . 4) ( ^) <* 


£=0 


(I.B.22) 




(r.x) = 3* T eS 2 ) (r,x-e) 


dq"U) 


k 


£>=0 


D h 


+ F 1 T 0^ 2) (r,x-^)(— | r -)d^ 


dq"U) 




(I.B.23) 
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The mixed cases follow in an analogous manner. Since the 
integrations are performed on the variable |, manipulations 
with respect to the radial coordinate r can be performed 
on the integrating kernels and the wall temperatures and 
heat fluxes together with the mixed-mean temperatures are 
immediately obvious. 

While the evaluation of the integral expressions for 
the arbitrary boundary conditions is considerably more 
difficult than the linear superposition of the previous 
examples, if we are provided with the fundamental solutions 
in sufficient detail, it can certainly be carried out. 

I.C. Objective 

The objective herein is to obtain the fundamental 
solutions for the situation where the fluid in the annulus 
has constant properties, is in fully established laminar 
flow, and the effects of internal energy generation and 
axial conduction are negligible. 

The fundamental solutions are sought in two ways; first, 
by solving (i.B.l) subject to each of the fundamental bound- 
ary conditions, and second, for the fundamental solutions 
of the second kind, through experimental measurement. 

In order to provide a complete description of laminar 
flow heat transfer in annuli these solutions are sought for 
a sufficient number of specific values of the radius ratio, 
r* , to encompass most practical sizes. 

I.D. Summary 

The methods described in sections II. A through II. H 
provide analytical expressions for the fundamental solutions 
and the numerical values necessary to evaluate these ex- 
pressions . 


14 



The values of the related non-dimensional quantities 
are presented in section II. I as functions of the axial 
distance parameter 

- x/D h 
x “ RePr ' 

for various values of the ratio of the radii, r * To 
facilitate interpolation on r* , the values of all the 
fundamental solutions obtained by McCuen 11 for the parallel 

plane channel are presented. 

A comparison of the fundamental solutions of the second 
kind with experimental values obtained as described in 
section III o A is presented In section III.B. 

Numerical examples of the use of the fundamental 
solutions in problems of superposition are considered in 
Appendix G. 

I,E. Conclusions 

The values for the fundamental solutions which are pre- 
sented provide a method for solving laminar flow heat transfer 
problems in annuli which is, essentially, complete. That is, 
values are provided for all x for enough values of r 
to enable the consideration of virtually any annular geometry. 

There are two effects which influence the departure of 
the behavior of the solutions from the simple, fully established 
(constant film coefficient) situation. These are. 

1. The effect of axial variation of boundary conditions; 

2. The effect of heating asymmetry, i.e., the influence 
of the boundary condition at the opposite wall. 

The axial effects will be significant near any discontinuities 
of boundary conditions or if the continuous variation is 
rapid. The asymmetric effects will predominate in long 
tubes and will be larger for a greater degree of asymmetry. 
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For cases of simple superposition like the examples on 
pages 10 to 12, probable criteria for the influence of these 
two effects can be presented as follows. 

For x < 0.01, the axial effects will predominate, 
for x > 0.1, the asymmetric effects will predominate, and 
for 0.01 < x < 0.1, the effects will be comparable. 

The good agreement obtained between the predictions of 
the theory and the measured values indicates the validity 
of the assumptions for fluids of Prandtl numbers near unity. 
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II THE FUNDAMENTAL SOLUTIONS 


II. A, Previous Studies 

Laminar flow heat transfer for the cases which represent 
the limiting forms of the annulus (i.e,, the circular tube 
and the parallel plane channel) has received considerable 
attention. 

Heat transfer to a fluid in laminar flow in a circular^ 
tube with constant wall temperature is the classical Graetz 
problem. The original work by Graetz and the extension of 
it by Sellars, Tribus and Klein 19 comprise a relatively 
complete solution corresponding to the limiting geometry for 
the fundamental solution of the third kind. Sellars et al 
also work out the case where the heat flux at the wall is 
constant by an inversion method. A direct solution for the 
fundamental solution of the second kind was obtained by 
Siegel, Sparrow and Hallman 20 . The fundamental cases, 
one and four, have no counterpart for flow in a circular tube. 

For laminar flow in a parallel plane channel all of 
fundamental cases have significance although they are symmetric 

fl( k ) = ) . All of the cases have received some 

' * , ii oo ' 

analytical attention. Perhaps the greatest effort has been 
for the case of constant and equal temperatures at the duct 
walls (sometimes called the Purday 1 ^ problem). More recently 
Cess and Schaffer have published solutions for cases one and 
two (Refs. 2 and 3). Case three has been considered by 
Schenke 18 . McCuen 11 has reduced all of the available solu- 
tions to a form consistent with the definitions given in 
section I.B and has completed the solutions of the third and 
fourth kind. 

The principal technique employed in obtaining these 
cases is that which Graetz used, solution of the energy 
equation by separation of variables. This method expresses 
the solution for the temperature field as an infinite series. 
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the terms of which are products of an exponential involving 
the axial coordinate x, the eigenfunctions of the resulting 
Sturm-Liouville equation, and a coefficient resulting from 
an orthogonal expansion of the initial condition. This 
approach has two main disadvantages: 

1. Near the step change in boundary conditions (x -» 0) 
the infinite series converges rather slowly 
necessitating a large number of terms. 


2. The eigenfunctions corresponding to the higher modes 
of the Sturm-Liouville equation are difficult to 
calculate directly. 

Q 

To avoid the series solution Leveque used the fact that very 
near the step change in boundary conditions the temperature 
signal has affected the fluid only in the immediate vicinity 
of the wall. He was able to obtain a similarity solution for 
the temperature field which becomes exact as x approaches 
zero. Sellars et al have applied the WKBJ approximation to 
the Sturm-Liouville system to provide an asymptotic form for 
the higher modes and to a large extent have circumvented the 
difficulty associated with that calculation. 

For true annular flow very few of the problems have 
received much attention. Jakob and Rees^ obtained the 
temperature distribution valid as x tends to infinity for 
the fundamental solution of the second kind. Singh‘S pro- 
posed a method for solution for a more general form of 
Eq. I.B.l which includes the effects of axial conduction and 
viscous generation. This method involves the expression of 
the eigenfunctions as a series of Bessel functions satisfying 
the boundary conditions. Murakawa 12 * has presented an 

integral equation formulation for the solution of the first 
kind and a series solution approach to the same problem 
including the more difficult case where arbitrary peripheral 
variations are allowed. Murakawa does not carry his solutions 
to the point of numerical computation. 
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II. B The Hydrodynamic Problem 

The velocity profile for fully established laminar flow 
in a concentric annulus with constant fluid properties was 

O 

first presented by Lamb . If we let 

r = r/r Q U = U/U m and r* = r ± /r 0 

where U , the bulk mean velocity is given by 
m 

U = , / r ° U(r)r dr (II.B.l) 

r 2 - r? r, 
o i i 

Then we may write 

U = | (1 - r 2 + B in r ) (II.B.2) 


where 


B = (r* 2 - l)/4n r* and M = 1 + r* - B 


The friction factor f is defined by 

t2 


T.r. + t r , 

f = — / ip J. 


r i + r o 


m 


(II. B. 3) 


where t, and t are the stresses exerted by the boundary 
1 o 

walls on the fluid at the inner and the outer surfaces respec- 
tively. The friction factor is related to the axial pressure 
gradient by 


dP 4f U 2 
„ m 

' dl = p "2" 


(II.B.4) 
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From the definition of the shear stress in laminar flow and 
Eq. (II.B.2) we can obtain 


_ 16 (l - r*) 2 
Re 


(II. B. 5 ) 


When r* = 0, the friction factor becomes 16/Re, the value 
for the circular tube. At r* = 1 successive application of 
L' Hospital's rule yields f = 24/Re, the friction factor for 
flow between parallel planes. Figure II.B.l presents f*Re 
as a function of the radius ratio r*. 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.^ 1.0 

r* = r./r 
i o 

FIG. II.B.l Friction factor for fully-developed 
laminar flow in a concentric annulus 


20 




II. C. The Energy Equation 


Let us introduce the non-dimensional variables 

- X / D h 

r = r/r Q , r* = iyr 0 x - RePr 


where 


D 


h 


- - 




into Eq. (I.B.l). The velocity distribution, U(r), is 
given by (II. B. 2). By virtue of the defining relationships 

we have 


d 2 d ( . k) 1 39 

J -f — u 

Sr 2 r 2r 


(1 - r 2 + B in r) 30^ 
2M (l - r*) 3x 


(k) 


(n.c.i) 


where, as before 


B = 


*2 n 
r* - 1 

in r* 


and M = 1 + r* - B 


This is the differential equation which the fundamental 
solutions must satisfy. 
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II. D., .Asymptotic (large x) Values of the Fundamental 

Solutions, Homogenization of the Boundary Conditions 

We will seek solutions to Eq. (iI.C.l) by the method of 
separation of variables. This method results in expression 
of the solution as an infinite series. In order to satisfy 
the initial condition we will need to express a function 
(representing the temperature distribution at x = 0) as an 
Infinite series of the orthogonal eigenfunctions. The resul- 
ting Sturm-Liouville equation will generate such a set only 
if the dependent variable satisfies general homogeneous 
boundary conditions. We note, however, that the boundary 
conditions specified for the fundamental solutions (l.B.4, 
I.B.6, I.B.8, I.B.10) are Inhomogeneous . We shall resort 

to the method of Siegel, Sparrow, and Hallman 20 . If we let 


_(k) 

e j 


= 0 


(k) 

j 


k) 

- ) 

j 'id 


(II.D.1) 


(k) 

where 0 is the "fully developed" temperature profile, 

that is, the asymptotic value of 0^ k ^ as x tends to 
infinity, then the will satisfy homogeneous boundary 

conditions . 


1* The fundamental solutions of the first kind 

We have defined ejj 1 ) as solutions to Eq. ( II.C.l) 
satisfying the boundary conditions: 

= 1 on wall j 

= 0 on the opposite wall, for all x > 0 
and = 0 for all r at x < 0 . 
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As x approaches infinity we see that, because the two walls 
are at different temperatures, whatever heat crosses wall i 
must cross wall o and hence /dx -» 0. Then Eq. (H.C.l) 

reduces to 


a 2 eb) 1 se' 1 ’ 

J 4- — J 

dr r dr 


0 


(II. D. 2) 


so that 


e 


(i) 

j 



In r + 


(II.D.3) 


Evaluating G ] _ and H 1 by the boundary conditions, we obtain 


fl (l)^ = lH?. (II. D. 4) 

1 jfd " in r* 


and 


e M\ = 1 _ (II. D. 5) 

o j fd in r * 

Now if we let = 9^ - 9^ ) fd we see that is 

a solution to equation (II. 0.1 ) and has homogeneous boundary 

conditions at r = r* and 1. That is, ~$ ± l (r,x) satisfies 
Eq. (II.C.1) and the boundary conditions 

^^(l,x) = 0 


^(r^x) = 0 


0pHr,o) = 


in r 
in r* 
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Similarly ( r * x ) satisfies Eq. (iI.C.l) and the 

boundary conditions 

^ l} (l,x) = 0 

^ X) (r*,x) = 0 

The asymptotic values of the can he obtained 

from Eqs. (II.D. 4 ) and (lI.D.5). 



*(l)) _ _ 2 (1 - r*) 

ii ; fd * „ * 

r * in r* 

(II. D. 6) 


®(D) _ a(i - r ») 

oi 'fd in r* 

(II. D. 7 ) 


*(!)) - . g(l - r*) 

00 ^fd in r* 

(II. D. 8) 


<1.(1)) _ 2(1 - r*) 

io 'fd r* in r* 

(II. D. 9 ) 

We may determine the asymptotic values of 0 ^}) and 

mi 

as follows. 


0 (k) 

mj 

= 7 2 J 1 U r dr 

[l - (r*) 2 ]r* J 

(II.D.10) 


Substituting, successively, Eqs. (II.D. 4 ) and(lI.D.5) into 
the above and carrying out the integration yields 



(II.D. 11 ) 


in r* 


q( 1 ) \ _ 4 

mi 'fd " M In r* 
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(II.D.12) 


e (l)) 4 

mo 'fd “ 1 M in r* 


B 3 7 ( r *}2 (r»)_ 

¥ " IF IF } + 4B 


Recalling the definition of Nusselt number, this can be 
related to the fundamental solutions as follows. 


Nu 


oo _ 

« ' k 


4> 


00 


- V 1 


(II.D.13) 


Then we can obtain 


Nu ^ ■*" ^ ) 
iNU ii ' fd 


il = Nu ‘° ))fd (II - Da4) 

L 2 " 8 8^ r ‘ + 2B 


Nu ^ O ) = 

w oo 'fd 


M(1 - r*) 


q + £(r*)' 


B 

2 


(rM\ 

2B J 




(II. D. 15) 


2. The fundamental solutions of the second kind 

We have defined as solutions to Eq. (H.C.l) 


J - = 1 on wall j. 


j 

satisfying the boundary conditions: 

( 2 ) 

= 5 *~ 

= 0 on the opposite wall, for all x 
and = 0 f° r a11 r at x < 0* 

As x approaches infinity d0j^/dx becomes constant and 
is equal to bO^j^/bx. Then Eq. (H.C.l) becomes 


3F 2 


i be 


+ — _ 
r dr 


( 2 ) 


tbe^ 


mj 


dx 

25 


(l - r 2 +B In r) 

— 5— (II.D.16) 

2M (l - r*) 



Equation (II.D.16) may be Integrated directly to yield 


3 'fd 


i^V ] 

m.1 

1 

"u - b) _ 2 

1 bx 1 

2M( 1 - r*) 2 

4 ’ 16 


Br 

+ 


+ Gg in r + 


(n.D.17) 


as the general form of the asymptotic value of the fundamental 

solution of the second kind. Prom the concept of 0^^ as 

a non-dimensional temperature we can evaluate d0'?V3x 
. mj ' 

from energy considerations. We have 


be 


(2) 

mi 


4r* 


dx 


(l + r*) 


(II.D.18) 


and 


c>0 


( 2 ) 


mo 


6x 


(1 + r*) 


(II.D.19) 


from which we can obtain 


,(2) _ 4r* 


mi 


(1 + r*) 


(II.D.20) 


and 


,( 2 ) 

mo 


(l + r*) 


X 


(II.D.21) 
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Using Eqs. (II.D.17) and (II.D.18) together with the appro- 
priate boundary conditions 


i ; fd 


(l + v*) M(l - r* )‘ 


U - Sl (r 2 - In ¥) 
2 


- f + | F 2 In r - + h r - *) - - kg 


( 2 ) 


+ e- 


(II.D.22) 


Here 0^2) is still unknown hut it can he evaluated hy 
noticing that 0^ - 0^' is a function of r* only as 
x gets large. Using Eq. (II.D.10) we can obtain 


9 (?> 


ei?) } 


ml fd (1 + r* ) M 2 (l - r* )‘ 


(-1 

73 

- 21 B 

+ U B 2 - 

1 

l 

.48 

18 

16 

2B_ 


25 2 5 r 2 

W “ 9 B ' IF B 


( r *) 2 - 


19 „ 11 

35 B - m 


/ „ \ 4 , 11 / *\6 

(r*) + Tfg-Cr*) 

(II.D. 23 ) 


Using Eqs. ( II.D.17) and (II.D.19) we can obtain a similar 


expression for 9 ^ ') 


0 (2)) = ± 

0 fd (l + r*) M(l - r*) 

2 


h - si p 2 - 1) 

2 


ix? - 1) 


8 


+ in F 


bull [ (r*) 2 
2 


- B] in r ) + 9 


( 2 ) 

OO 


(II.D.24) 
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,( 2 ) 


As before can be evaluated from (II.D.24) and (II.D.10) 


e (2) . e (2)l = 1 

00 mo fd (1 + ~) M 2 (l - r*) 2 


— - — B + — B 2 


48 36 


16 


+ 


11 , 2B 11 „2 

¥8 + 9~ “ T6 B 


(r*)‘ 


25 31 -□ 

W ‘ 18 B 


(r*) 4 - £f(r*) ( 


+ lull 

+ 2B 


8 


(II.D.25) 


Now If we let = 9^^ - we see that 

J J J fd j 

Is a solution to Eq. (II.C.l) and has homogeneous boundary 

conditions at r = r* and 1. Explicitly, 15^ 2 )(r,x) 

satisfies Eq. (II.C.l) and the boundary conditions 


W 


( 2 ) 


— — — =0 at r = r* and 1, for all x, and 
dr 


S^ 2 )(r,o) = 


(l + r* ) M( 1 - r*)‘ 


(r 2 - to 7) 

2 


. f. + | F 2 tn - . (r>) g M + (1 „t g to r*) _ ^ 


- (*<?> - 8 <?>> 


ii mi 


/ 2 ^ 

and ^ '(r,x) satisfies Eq. (II.C.l) and the boundary 


o 

conditions 


sr(2) 


be 


dr 


= 0 at r = r* and 1, for all x and 
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tf^(r,o) = ± 

° (1 + r*) M(l - r*) 


1 - B) t—2 


[r 4 - 1) 


+ 5 r 2 In r - [(r*) 2 - B] in r 


' oo mo ' f d 


Using Eqs . (II.D.22) and (II.D.23) we can obtain the 
asymptotic expression for 


9 (2) . e (2)) 

oi mi ' fd 


(1 + r*) M 2 (l. - r*) 


5 11 -d i 11 -q2 1 i 

8 ‘ T B + TS B " 25 j 




+ B " ^8 1 (r#)4 " M (r * )6 


1 + in r*] 
2 
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and similarly from (II.D.24) and (II.D.25) we get 


»(2) 

io 


e (2 >) 

mo 'fd 


(1 + r*) IT(l - r*)‘ 


-1 + ^ 


_5 B 2 


48 36 16 


13 , 5 -D . 11 -o2 
Iff + IS B + T6 B 


(r*)‘ 


+ 


17 31 

W ” 18 


B 


( r *) 4 + *§ ( r *) 6 - ^ 


8 


(r*) in r* 
2 


(II.D.27) 


Because of the specified boundary conditions we see 

( 2 ) f 2 ) 

that Nu oj/ = Nu^ ' = 0j and 


Nu 


( 2 ) 

ii 


fl (fT : e (?) 

ii mi 


(II.D.28) 


Nu 


(2) » 1 
00 Q { 2 ) _ e ( 2 ) 


00 mo 


(II.D.29) 


3 • The fundamental solutions of the third kind 

We have defined 0^3 ) as solutions to Eq. (iI.C.l) 
satisfying the boundary conditions 

0p) = 1 on wall J 

= -D^ d0^V^ n = 0 on the opposite wall, for all 
x > 0 
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and 


,(3) _ 


= 0 


for all r when x < 0. 


As x approaches infinity the temperature of the fluid 
tends toward a constant equal to that of wall j. Hence 


q(3)> = 9^3)) = 0^)) = 0^3)) 

0 i 'fd y ii 'fd oi 'fd mi 'fd 


= 1 


(II.D.30) 


<$( 3 ) ■) = <x>(3) ) = 0 

®ii ' fd oi ' fd 


(II.D.31) 


and 


<4 3) > 


fd 


0 ( 3 )) = g(3)\ _ g(3)> = 1 

0 oo 'fd y io 'fd y mo ^fd 


(II.D.32) 


<j(3) ) = ^3) ) = 0 

00 'fd io 'fd 


(II.D.33) 


Now if we let 0j 3 ^ = ” ®j^^fd we see 

satisfies Eq.(lI.C.l) and has homogeneous boundary conditions 
at r = r* and 1. Explicitly 0^ 3 ^(r,x) satisfies 
Eq.(lI.C.l) and the boundary conditions 


0(3) ( r * jX ) = 0 


<J>( 3 )(l,x) = 0 and 0^ ; (r,O) = - 


a^3 ) ( „ 


and 0^ 3 ^ (r>x) satisfies Eq. (iI.C.l) and the boundary 


conditions 
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e£ 3 ) U,x) = o 

■^(3)( r *^— ) _ o and "0^ ^(F,0) = -1 

The limiting Nusselt numbers reduce to an Indeterminate 
form and cannot be evaluated from these considerations. 

4 . The fundamental solutions of the fourth kind 

We have defined the 0j^ as solutions to Eq. (iI.C.l) 
satisfying the boundary conditions 

(4) Se $ 4) 

" D h ”3n — = 1 on wall j 

0^^ = 0 on the opposite wall for all x > 0 

(u) - - 

and 0j ' = 0 for all r when x ^ 0 

As x approaches infinity the heat passing through 
wall i must also pass wall o and ^0j^V^x -* 0. Then 
0j^^ will satisfy equation (II. D. 2) and 

0 j 4)) fd = g 4 ^ * + H 4 (II.D.34) 


By application of the appropriate boundary conditions one 
obtains 


Q W) 

i 'fd 


2(1 - r*) 


in r 


and 


e^h = I 

0 fd 2(1 - r*) 


in (£ ) 


(II.D.35) 


(II.D. 36 ) 
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Now if we let 0^ = 0^ - 0^ ) - , we see that 

( ii \ .3 V J J ia 

0^ ' satisfies Eq. (II.C.1) and has homogeneous boundary 

conditions at r = r* and 1. That is ej 4 '(r,x) satisfies 

Eq. (II.C.1) and the boundary conditions 

$( 4 )(r*,x) = 0 

■0^ 4 )(l,x) = 0, 0^(r,O) = 2 (r r z~ r ~ r j ^ r 

and 0^ 4 )(r,x) satisfies Eq. (iI.C.l) and the boundary 
conditions 

^ 4) (l,x) = 0 

^o^(r*/x) = 0 

Then as before 


and 0^(r,O) = " 2(1 - r*) ^r*^ 


4 ) \ _ -r* l n r* 

ii ; fd 2(1 - r*j 


(II.D.37) 




-r 


* 


(11.D.38) 


£ ,<*>) 
r\r\ * 


00 ' fd 


Jt n r* 
2(1 - r*) 


(II.D. 39 ) 


and 


*<*)) 

*io 'fd 


(II.D.40) 


Further 

fl (4U 2r* 

0 mi 'fd ~ MU - r*; 


h + ^ 


B (r*) 
“ T ‘ "4B 


(II.D.41) 
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(II.D.43) 

(II.D.44) 


The large x values of the various non-dimensional 
quantities are given, as functions of r*, In tables II.D.l 
to II. D. 3 and In Figs. II.D.l to II.D.10. 

The nature of many of the algebraic expressions Is such 
that when r* -» 1, they present an Indeterminate form. 

At r* = 1, this can be avoided either by using L 1 Hospital's 
rule or by transforming Eq.(l.B.l) to rectangular coordinates 
and generating the solutions to the transformed equation. 

For r* near unity, however, the evaluation of many of the 
preceding forms involves the differences of nearly equal 
quantities with a resulting loss of accuracy. This is 
reflected in the tables by the presentation of only three 
significant figures for r* = 0.9 (r* = 0.8 and 0.9 for the 
second kind). 

Examples of the fully-developed temperature profiles 
are given in Figs. II.D.ll to II.D.14. 
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TABLE II.D.l 


ASYMPTOTIC (LARGE Y) VALUES OP THE FUNDAMENTAL SOLUTIONS OP THE FIRST KIND 


r* 

*(U. _*(1) 

11 T lo 

*(D. .*(1) 

00 ol 

0 U) 

9 ml 

a (4 

mo 

““io 5 

00 01 

0 

00 

0 

0 

1 


2.66667 

0.02 

25.0510 

0.501020 

0.169969 

0.830031 

30.1808 

2.94772 

0.04 

14.9120 

0.596482 

0.198875 

0.001125 

18. 6139 

2.99928 

0.05 

12.6847 

O.634236 

0.210090 

0.789910 

16.0584 

3.01887 

0.06 

11.1371 

0.668228 

0. 220072 

0.779928 

14.3042 

3.03640 

0.08 

9.10607 

0.728486 

0.237760 

0.762240 

11.9465 

3.06395 

0.10 

J. B 1730 

0.781730 

0.252556 

0.747444 

10.4587 

3.09528 

0.20 

4.97068 

0.994136 

0.309374 

0.690626 

7.19736 

3.21338 

0.25 

3.32809 

1.08202 

0.33H97 

0,668803 

6.47139 

3.26700 

0.30 

3.87606 

1.16282 

0.350455 

0.382846 

0.649545 

5.96629 

3. 31803 

0.40 

3.27407 

1.30963 

1.44270 

0.617155 

5.30510 

3.42077 

0.50 

2.88539 

0.409987 

0.590013 

4.89038 

3.51888 

0.60 

2.61015 

1.56609 

0.432800 

0.567200 

4.60183 

3.61851 

0.70 

2.40315 

1.68220 

0.452770 

0.547230 

4.38978 

3.71536 

0.80 

2.24071 

1.79257 

0.470300 

0.529700 

4,23014 

3.811540 

0.90 

2.11 

1.90 

0.466 

0.514 

4.11 

3,9i 

1.0 

2.00000 

2 . 00000 

0.500000 

0.500000 

4.00000 

4 . 00000 


ASYMPTOTIC (LARGE X) VALUES OP THE FUNDAMENTAL SOLUTIONS OP THE SECOND KIND 


f a(2) fl (2h 

6 ml * 


0. 0561442 

0.0627609 

0.0742498 

0.0839928 

0.117662 

0.128974 

0.138099 

0.151899 

0.161785 

0.169154 

0.174759 

0.179 

0.183 

0.185714 


in (2) e { 2 ) ) 

' 9 oi " 9 ml J 


-0.00255890 

-0.00496486 

- 0.00612638 

- 0.00726861 

-0.00948471 

-0.0116214 

-0.0212779 

-0.0255187 

-0.0295091 

-0.0365981 

-0.0427568 

-0.0461504 

-0.0529590 

-0.0570 

-0,0609 

-♦0.0642857 




0.229167 

0.211227 

0.209292 

0,208682 

0.206193 

0.207439 

0.206859 

0.204809 

0.203884 

0.202922 

0.200836 

0.198549 

0.196108 

0.193554 

0.191 

0.188 

0.185714 


-0.145833 

-0.127945 

-0. 124122 

-0.122567 

-0.121143 

-0.11B559 

-0.116214 

-0.106390 

-0. 102207 

-0.0983635 

-0,0914955 

-0.0855135 

-0.0802490 

-0.0755920 

-0.0717 

-O.0677 

-0.0642857 


32.7051 

20.5092 

17. 8113 

15.9335 

13-4681 

11.9058 

8.49692 

7.75348 

7.24116 

6.58331 

6.18102 

5.91176 

5.72215 

5.58 

5.47 

5.38461 


4.36364 

4.73424 

4.77802 

4.79198 

4.80322 

4.82070 

4.83421 

4.88259 

4.90475 

4,92801 

4.97918 

5.03655 

5.09922 

5.16653 

5.24 

5.31 

5.38461 


ASYMPTOTIC (LARGE Y) VALUES OF THE FUNDAMENTAL SOLUTIONS OP THE FOURTH KIND 



e (4) 



r 

9 11 

oi 

00 

0 

0 

0 

ce 

0.02 

0.039919 

-0.02 

1.99593 

0.04 

0.067060 

-0.04 

1.67650 

0.05 

0.078835 

-0.05 

1.57670 

0.06 

0.089670 

-0.06 

1.49650 

0.08 

0.109017 

-0.08 

1,37271 

0.10 

0.127921 

-0.10 

1.27921 

0.20 

0.201180 

-0.20 

1.00590 

0.25 

0.231049 

-0.25 

0.924196 

0.30 

0.257994 

-0.30 

0.859979 

o.4o 

0.305430 

-0.40 

0.763576 

0.50 

0.346574 

-0.50 

0.693147 

0.60 

0.383119 

-0.60 

0.638532 

0.70 

0.4 16209 

-0.70 

0.594584 

0.80 

0.446287 

-0.60 

0.557859 

0.90 

0.474 

-0.90 

0.527 

1.0 

0 . 500000 

-1.00 

0.500000 


9 ^ Nu^-Nu^ 

mo ii io 


- 50.0000 
- 25.00000 
- 20 . 00000 
-16,66667 
- 12.50000 
- 10.00000 
- 5.00000 
-4,00000 
-3.33333 
- 2.50000 
- 2.00000 
-1.66667 
-1.42857 
-] .25000 
- 1.11 
- 1.00000 


0 

0.006785 

0.013337 

0.016563 

0.019760 

0.026110 

O.O32307 

0.062240 

0.076523 

0.090385 

0.116933 

0.142091 

0.165814 

0.168407 

0.209889 

0.230 

0.250000 


I.65669 

1.34309 

1.24545 

1.16716 

1 . 04663 

0.956141 

0.694700 

0.618104 

0.558695 

0.471244 

0.408966 

0.362175 

0.325431 

0.295499 

0.272 

0.250000 


30.1808 

10.6139 

16.0584 

14.3042 

11.9465 

10.4587 

7.19736 

6.47139 

5.96629 

5.30510 

4.89038 

4.60183 

4.38978 

4.23014 

4,11 

4.00000 


2.66667 

2.94772 

2.99928 

3.01887 

3.03640 

3.06395 

3.09528 

3.21338 

3 . 26700 
3.31803 
3 . 42077 
3. 5188B 
3.61851 
3.71536 
3.81154 
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ASYMPTOTIC (LARGE 5?) VALUES OF THE ASYMPTOTIC (l_ARGE x) VALUES OF THE 

FUNDAMENTAL SOLUTIONS OF THE FIRST KIND FUNDAMENTAL SOLUTIONS OF THE FIRST KIND 



FIGURE n. D. I FIGURE II D. 2 


















PTOTC (LARGE X) VALUES OF THE ASYMPTOTIC (LARGE X) VALUES OF THE 

lIMTAL SOLUTIONS OF THE SECOND KIND FLMDAMENTAL SOLUTIONS OF THE SECOND KIND 



FIGURE H. D. 6 



TOTIC (LARGE 5?) VALUES OF THE ASYMPTOTIC (LARGE X) VALLES OF THE 

NTAL SOLUTIONS OF THE FOURTH KIND FUNDAMENTAL SOLUTIONS OF ThE FOURTH KIND 
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FIGURE H. D. 7 FIGURE H. D. 8 






TOTIC (LARGE X) VALUES OF THE ASYMPTOTIC (LARGE X) VALUES OF THE 

NTAL SOLUTIONS OF T>€ FOURTH KIND FUNDAMENTAL SOLUTIONS OF THE FOURTH KIND 



© 


FIGURE H. D. 9 FIGURE E D. 10 
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Fully developed temperature profile for a 
fundamental solution of the first kind 

FIG. XI. D. 11 



Fully developed temperature profile for a 
fundamental solution of the second kind 

FIG. II.D.12 
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1.0 






Fully developed temperature profile for a 
fundamental solution of the third kind 

FIG. II.D .13 



Fully developed temperature profile for a 
fundamental solution of the fourth kind 

FIG. II.D.14 
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IX. e. Limi ting (small x) Leveque Approx i mations to the 
Fundamental Solutions 

We seek a solution of (I.B.l) valid near the step change 
in boundary conditions which accompanies the fundamental 
solutions. Following the scheme of Lev/eque we will replace 
the expression for the velocity by a linear profile having 
the same slope at the wall. Differentiating (II.B.2) we 
have 


dU 

■Sr 


r=r 4 


2U 


m 

r M 
o 



2r*) 


(II.E.1) 


and 


SU 

■3? 


2U 

__ m 

" r M 


r=r Q o 


(B - 2) 


(II.E.2) 


Accordingly let us use as the velocity near the inner wall 


U i = 


2U 

m 

~W 


(S _ 


2r*) ( r 


- r 1 


(II.E.3) 


and near the outer wall 

2U 

U Q = -jjp (B - 2) (r - 1) (II. E. 4 

Consider first the solutions where the step change in 
boundary conditions occurs at the inner wall. Substituting 
(II.E.3) Into (I.B.l) and transforming to the non-dimensional 
variables } r and x we have 
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(II. E. 5) 


d 2 e 

hr 2 


+ 1 = E (r - r* ) 

r dr dx 


where = N(B/r* - 2r*) 
and N = 1/2 M(l - r*) 2 

as the equation which the fundamental solutions must satisfy 
near the step change in boundary conditions. 

This equation, and the analagous one obtained by using 
(II.E.4) in (i.B.l) near the outer wall, will be valid only 
for x sufficiently close to the step change in boundary 
conditions (x = 0) so that in the radial distance for which 
the temperature deviates from the initial condition, 0=0, 
the velocity profile is closely approximated by (II. E. 3 ) 
or (II.E.4). 


1. The fundamental solutions of the first and third kind 

Since the temperature signal only affects the fluid 

in the immediate vicinity of the inner wall, the boundary 

condition on the outer wall is immaterial and a step change 

in wall temperature provides the small x approximations 

to the solutions of the first and third kinds. Further for 

the region close to the wall the effect of curvature is slight 
1 bd 

and the term — — can be neglected and Eq„ (II.E.5) becomes 
r 5r 


d 2 0 

3F 2 




(ll.E.6) 


We now seek a similarity solution of the form 
0 = f(p) where p = ^ ~^ r *^ 
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Substitution into (II.E.6) yields the ordinary differential 
equation 

f" + a p 2 f' =0 (II. E. 7) 

where a = E^/3 

The boundary conditions in terms of f(p) become 

f = 1 at p = 0 

f = 0 as p -» 00 

The solution to (II.E,7) is 

3 

f = G / e 3 ~ dp + (II. E. 8) 

Applying the boundary conditions we obtain 


G 1 « ap^ 

/ e- -r dp 

o 

and o 

_ 

H x = 1 - [0 X / e 3 dp] p _ Q 


We need to evaluate 



but 


00 3 i 

/ e~° da = (i) I = 0.89298 
o J 


(II. E. 9) 


(II.E.10) 
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Then 


G 


1 


-1.1198 



Now we can write 


4 k) (r jX ) 



1 - 1.1198 



o 


(II.E.11) 


(II.E.12) 


for k = 1 or 3. 

We are Interested In the value of . From the 

definition of $ we get 



2(1 


r* ) 


d0\ 

~bb'o= 0 


bo 


(II.E.13) 


Then from ( II.E.12) 


9 


(k) _ 


11 


= 2.2396(1 - r*) 


B 0 „ 

7* - 2r* 

|_18(l - r*)^Mj 


nl/3 


(x)- 1/3 


(II.E.14) 


for k =1 or 3* 

The variation of mean temperature can he deduced from 
energy considerations. We have 


^ fi (k) 

° th* _ $ (k) 

bx 1 + r* lj - 


(II.E.15) 
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from which 

M _ 


9 


ml 


l6r* \ 

fv - 2r* 

2.2396(1 - r*) 

/ % 2 

\l+r*j 

L18(1 - r*; MJ 


1/3 



(ii.e. 16 ) 


for k = 1 or 3. 

2. The fundamental solutions of the second and fourth kinds 

Again since the temperature change is not felt at the 
outer wall, a step change in $ provides the approximation 
to the solutions of the second and fourth kinds. 

Because the derivative at the wall must be independent 
of x , let us seek a similarity solution of the form 
e = £ n f(p), where p = (r - r*)/^, and £ = x. 

Equation (II.E.5) becomes 

f"(p) + n E ± x(3n-l) p 2 f , ( p ) _ n E ± x^3n-l) p f(p) = 0 
If n is chosen as l/3.> x vanishes and we have 


f n (p) + a p 2 f'(p) - a p f(p) = 0 


(II.E.17) 


where a = E^/3. 


The solution to Eq. 


f(p) = G 2 P I 


(II.E.17) is 

_ 

3 

2—X — dp + H 2 p 


(II.E.18) 
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Then 



9 = 0 when x -* 0 ( p -» 00 ) , 

and 

M ) = _ 1 

dr r=r* 2( 1 - r*) 

We see that 

1 

^ ) =x 3 f'(p)] (%)= f'(p)) (II.E.20) 

dr r=r* / p=0 ’dr / p=0 

From (II.E.19) we obtain 

3 

r - Ml -I 

f'(p) = H - G ? /ape 3 dp 

/ P=0 i d -I p=0 

The integral on the right vanishes which can be seen by ■ 
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expanding the exponential in a series and integrating term 
by term. Then 


H, 


2 “ 2(1 - r*) 


(II.E.21) 


The remaining boundary condition demands that 


0 = (r - r*) 


2TT rFJ - u 2 


3 

_ afil 

M / a pe 3 dp) p=oo 


(II.E.22) 


Since (IF - r*) ^ 0, the bracketed term must vanish. We have 
already seen that the integral in (II.E.22) vanishes at p = 0, 
so we may write 


(a / pe 


3 

ap J 

3 


dp) 

•^ / p=oo 


a 0 J pe 


3 

ap^ 

3 


dp = 


1 2 

££ ( 2 ) , 

2 '3' * 


From this 


q = -0.36925 

2 (1 - r*)^3j a/3 ' 


and 


e[ k) (r,x) = (x) 3 


H 

0.36925 


2(1 - r*) 


(II.E.23) 


L(1 - r *)-^ a/3 1 


ap^ 

(e 3 + p / ap e 


ap^ 

3 dp) 


(II.E.24) 


for k = 2 or 4. 

We are interested in the value at the heated wall. 
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6 


(k) _ 0.36925 
II (l - v*) 



18(1 - r*) 2 MJ 


1 



(II.E.25) 


for k = 2 or 4. 

If the step change in boundary conditions occurs at the 
outer wall, U is given by Eq. (II.E.4) and Eq. (II.E.6) 
becomes 


dr 2 


= E o< 7 - 1>|§ 


where E q = N(B - 2). 


1 

If we now define p = (r - l)/x Eqs. (II.E.7) and (II.E.17) 
result where a = E q / 3 . Then, with the exception of the 
expression for the mean temperature In cases 1 and 3, the 
only difference is in the factor a , and we have 


^ = - 2.2396(1 - r*) 


B - 2 


18(1 - r*)^M 


_1 

3 


x ' 3 


(II.E.26) 


and 


9 


(k) _ 


mo 


6 


\! + r*, 


2.2396(1 - r *) 


B - 2 


18(1 - 


r*) 2 M 


- 1 P 

T £ 

3 - 3 

X U 


(11.E.27) 


for k = 1 or 3. 


Further 

q { k) = 0.36925 

°° " " (1 - r*) 

for k = 2 or 4. 


B - 2 


18 (l - r*)^M 


1 1 
3 r 3 


x 0 (II.E.28) 
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All of these equations have the form of a constant 
multiplied by a power of x. Table II.E.l presents these 
constants for various values of r*. They are given graphi- 
cally in Figs. II.E.l and II. E. 2. The symbols used are as 
follows . 


$ 


0 


00 

n 

oo 

mj 


- 1 



2 

C^) x ^ for j = i or o, k = 1 and 3 
mj 


and 




1 

3 


for J = i or o, k = 2 and 4. 
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II. F. Solution of the Energy Equation 

The will be solutions to Eq. (iI.C.l) satisfying 

the homogeneous boundary conditions as outlined in 
section II. D. Omitting the sub- and super-scripts for brevity 
let us assume that 6 = R(r)X(x), where R(r) and X(x) are 
functions only of r and x respectively. Then by the method 
of separation of variables we obtain the two ordinary differen 
tial equations 


X f + A 2 x = 0 


(iI.F.l) 


and 


where 


R" + i R« + A 2 N(1 - r 2 + B in r)R = 0 (II. F. 2) 
r 


N = 1/2 M(l - r*) 2 
Equation (iI.F.l) yields 

, 2 - 

X(x) = Ce A x (II. F. 3) 

where C is an arbitrary constant. Equation (II.F.2) Is 
of the Sturm-Liouville type* and will, in general, give rise 
to an infinite set of eigenvalues A , n = 0,1,2,..., and 
corresponding eigenfunctions R , n = 0,1,2,.... It can be 
shown that these functions R ($•) form an orthogonal set 
with respect to the weight function 

w(r) = N(r - r 3 + Br in r) (II. F. 4) 

that is 1 _ ___ __ 

/ w(r) R (r) R (r)dr = 0 if n / m (II. F. 5) 

r * n m 

A comprehensive discussion of Sturm-Liouville theory may be 
found in Ref. (4). 



Further, this orthogonal set is complete so that we may ex- 
press an arbitrary function (having proper smoothness con- 
ditions) as an infinite series of the eigenfunctions as 
follows: 

F(r) = Z C R (r) (1I.F.6) 

n=0 n n 

where 

1 _ _ _ _ 

/ w(r) F(r) R (r)dr 

C = r* (II. F. 7) 

n 1 _ o _ _ 

/ w(r) R^ (r)dr 
r * 

In view of this we can write the desired solution as 

0'(r,x) = S C R n (r)e ^n x (II.F.8) 

n=0 

with the A n and the C n so chosen as to satisfy the 
appropriate boundary and Initial conditions. We are now 
able to determine the boundary conditions in terms of 
Eq. (II. F. 2) and to write an explicit expression for (II. F. 7) 
for the various fundamental cases. 

1. The first kind 

The boundary conditions for this case require that for 
both and 9J 1 

R n> ( r .) = R U) (!) . o (II. P.9) 

In essence this is the equation of the eigenvalues. The 
values of the parameter In Eq. (II.F.2) must be chosen 

such that the solutions generated satisfy the condition 
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(II.F.9) above. We see then that (A = (A The 

' n'i v n'o 

coefficients, C , in Eq. (II.F.8) can he obtained from 
(II.F.7) with F(r) = 0^ ^(r,0). With some manipulation of 
the differential equation (II.F.2) we can obtain* 



(II.F.11) 
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(II.F.14) 


3 . The third kind 

Here the boundary conditions oh the eigenfunctions are 
not symmetric so we must consider separately the cases where 
the temperature step is applied at the inner or the outer 
walls. 

For we have 

K (r*) - *a w , - o (II. P.15) 

n a? r_1 


and 



(II.F.16) 
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we have 


For 6^3) 
o 


SR i 
n 

l Sr I — 


and 


(C ) 
v n'o 


(3) _ 


= R (1) = 0 
n 


K\ 

1 7=1 


n 


R (r*) r* — 
n Sr 


(II.F.17) 


dR \ 

dR » 

7 dR 1 


n + 

11 

n 


l SA / - * 1 

isA 

Ur i 

, 

n r=r* 

n 

r=l 


Note that here (A )^) ^ ) ^ ) 

' n'i ' n'o 


4. The fourth kind 

—(L) 

For 0^ ' we have 


SR 


and 


(C )V 
v n'i 


(4) , 


a _ n )*. r * - R n (l) " 0 


R (r*) 


n 


R n (r*) r* — 
Sr 


(II.F.18) 


(II.F.19) 


r^n 

rsRH 


fSRH 

-] 

n 1 + 

n 


n 


, SA / - * 

ISA 


> Sr J 


n r=r* 

n 

r=l 


(II.F.20) 


since ( II.F.19) is identical with (II.F.17) we will have 

(a )$ 4) = (a ) (3) 

' n'i ' n'o 
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we have 


For 9 




dR 

R Jr*) =-=?)_ = 0 

dr r=l 


(II.F.21) 


and 


(C ) 
v n'o 


(4) _ 


R (1) 

n' ' 



r a 

i aR \ 


f dR \ 

A 

R n ( 1 ) — 

—a + r* 



n 

n dF 1 

idA - , 

IdA 

1 dr - . 



n r=l 

n 

r=r 


(II.F.22) 


Note that 


(A ) (4 > = (A )( 3) 
' n'o ' n'l 
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II. G. Behavior of the Solutions as Lambda Becomes Large. 

We would like to examine the behavior of the solutions 
to (II.F.2) as A becomes large. Following the method 

r IQ 

of Sellars , Trlbus and Klein , let us develop the so- 
called WKBJ*approximation to the solution of ( II.F.2). If 
we let 


we find that g(r) must be a solution to 

g"(r) + (g'(r)) 2 + &llll + a 2 n( 1 - r 2 + B In r) = 0 

r n 

(II.G.1) 

An asymptotic solution is sought in the form 


g = Ag Q + g 1 + A -1 g 2 + ... (II. G. 2) 

Substituting (II.G.2) into (iI.G.l) and observing that to 
satisfy the equation for all values of A n , the coeffici- 
ent of each power of A n must vanish, we can obtain 


g^ = ± i VN(1 - r 2 + B Tn r) 

and 



(II.G.3) 


(II.G.4) 


Since A^ is large, the remaining terms in ( II.G.2) are 
neglected. From (II.G.3) and (II.G.4) we may write 


After Wentzel, Kramers, Brillouln, and Jefferies. 
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G e 


R n (r) 


+ i A 


/ "y^Cl-r 2 + B in r J dr 
r 


+ H 


-i A 


e 


/•y^l-? 2 + B * n r> ' 



[n( 1-r 2 + B In r)] 1/A (II. G. 5) 

This is the WKBJ approximation. Equation ( II.G.5) must be 
patched unto the regular solution of (II.F.2) near the 
boundary walls. Equation ( II.G.5) can be put into a some- 
what more useful form for our purposes which is 


R n (r > = 


G cos 
o 


r 2 + B in r) dr - 0 
r* 


N(l- 


_2 

r + B in r ) J 


(II. G. 6) 


Near the inner wall, r -» r* . Let r = r* + z and 


o p p Bz 

(l-r 2 +Binr)=l- r* - 2r*z -z + r* - 1 + pr 


Since z is small we neglect terms of higher than first 
order and 

(1 - r 2 + B in r) z - 2r*) as z -- » 0 

Equation ( II.F.2) becomes 


R" + - R 1 + A 2 E z R 
n “ n n i n 


0 


where 


E ± = N(B/r* - 2r*) 
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Let 


2/3 


T] = A n w J z , then we have 


d 2 R 


n 


dR 


dT1 ^ A n 2// ^ ( v * + nA^' 5 ) dr] 


“ 273 ~ 


n 


+ E ± tj R n = 0 (H.G.7a) 


which, if is large, becomes 


d 2 R 


p 'H E — 0 ( II.G.7b) 

dr] 

Equation (lI.G.7b) is a form of Bessel's equation and has the 
solution 


R n * G 1 T)l/2 J l/3 ( -^ r|3/2) + H 1 1 1/2 T ' 3/2) 


(II. G. 8) 


We desire Eq.(lI.G.6) to correspond to Eq.(lI.G.8) as r -* r*. 
Observe that for this condition the integral appearing in 
(II.G.6) becomes 


o l/2 do - |V5; z 3/Z , andVT’-V^' 


Using this we have 


R 

n 


r-*r* 


G 


o 


cost 2/3 Ve; A n z 3 / 2 _ 0 ] 

V? 7 z 1/4 


(II.G.9) 
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For large t) , Eq. (II.G.8) can be written 


R = 
n 


1 

[ 


< G, cos 

5 A 1/6 

\ 1 

3 n - 



r , 3/2 !l 

3 n z '12 


+ cos 


rS ^l A z 3/2 . TL- 

3 n Z 12 


E 1 1>/4 z 1//4 (II.G.10) 


If G is taken as 
o 

a o - -^/Jl \ 1/6 (II.G.ID 

then Eqs . (II.G.9 ) and ( II.G.10) are Identical In form. 
Expanding the cosine of the difference of angles occuring 
In (II.G.9) and (II.G.10) one obtains as the equations for 
the coefficients 

G^ cos (5 tt/12) + H 1 cos(7t/12) = cos0 (II.G.12) 

G^ sln(57r/l2) + sln(7r/l2) = sin0 (II.G.13) 

Since 0 is, as yet, an unknown quantity, ( II.G.12) and 
(II.G.13) are insufficient to define completely G^ and 
H 1 , and we need to apply the appropriate boundary conditions. 
Let us reserve this to a later time and continue the general 
development . 

Near the outer wall, r -* 1. Let r = 1 - I and 
(1 - r 2 + p in r) -* 4(2 - p) 

Let £ = A £ , then for large A Eq.(lI.F.2) can be 
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reduced to 


d 2 R 


§ + E £ R =0 
2 o n 


where 


d C 


E = N( 2 - B) 
0 


(II.G.14) 


Equation ( II.G.14) has the solution 


*n - ^ V3 < ^ + ^ ^ 


(II.G.15) 


Here we desire (II.G.15) to correspond to (II.G.6) as r -» 1, 
For this condition the integral appearing in ( II.G.6) becomes 


r 2 + B in r) dr - 

r * 



where 


7 = 




— IT 1 _ 

+ B £n r)dr 


(II.G.16) 


Using this we have 

Vr 7 co s [ >^ n 4 3//2 - y + <t>] 

e lA ^lA 
V 3 n o 

For large £ Eq. ( II.G.15) may be written 


(II.G.17) 
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R = 
n 


.1 \ 1/6 

L3 n 


2\/F o/p 5 tt 
G rt cos [ ^s- 2 - A ^ 3/ " 12 ] 


n 


&nT o/p TT 
+ H 0 cos [ — 5 — — A„ - tt=t] 


n 


E 


lA ^lA (U.G.18) 


An appropriate choice of G 2 and H 2 would make (XI.G.17) 
and (II.G.18) Identical. As before, we obtain as the equa- 
tions of the coefficients 


G 2 cos(5tt/12) + H 2 cos(tt/12) = r* cos(y A r - 0)(II.G.19) 


G 2 sln(5ir/l2) + Hg sin(ir/l2) = r* sin(Y A n - <t>){ II.G.20) 


To complete this analysis we must apply the specific boundary 
conditions for the various fundamental cases. 


1 . The fundamental solutions of the first kind 


For this case the boundary conditions are given by 


the term 


(II. F. 9). In Eq. (II.G. 8 ), as q -* 0 
^i /2 j (1 ^ k ■q 3/2) does not vanish but approaches a 


constant. Then in order to satisfy the inner wall boundary 
condition H, must equal zero. From (II.G.12) and (II.G.13) 
we see that this demands that 


sin 0 cos(5tt/12) - cos 0 sln(5ir/l2) = 0 


or 

sin (0 - 5 tt/ 12 ) = 0 


(II.G.21) 
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( II. G. 21) will be satisfied if 

0 = if ± n7r * n = °; lj 2, ... ( II. G. 22) 

There Is no loss of generality If n Is taken to he zero. 
Further 


G 


1 


= 1 


(II.G.23) 


Then, near the Inner wall, we have 

R n = *„ l/3 , l/s J l/3 (^ A n z 3/ 2) 


(II.G.24) 


At the outer wall we must Insist that vanish, 
hence 


sin(Y A n - 0 - 5 tt/12) = 0 (II.G.25) 

(II.G.25) will be satisfied If 


A n = + §) f * n = 0, l, 2, ... (II.G.26) 

where 

7 = / ~^N(l - F 2 + B in r)'dr 

This is the equation for the eigenvalues as n becomes 
large. Further 


g 2 = (-i) n V^ 


(II.G.27) 
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Then for r close to unity we have 


R 


n - (-1)" -\P A n l/3 ? 1/2 J l/3 ( -j- 2 \ &'*) (II. 0.28) 


,3/2. 


We still need the expansion coefficients given by (II.F.10) 
and (II.F.11). Differentiating (II.G.24) we have 


dR 

— ) = K 

hr r-*r* 


,3/2. 


/o Jn /o ( — 5 A Z ) 1 

1/3 V3* 3 v . 2 A n ^, nZj4 (_i An 

3 


2Ve7 ,/p 

* z J,,M K z J/ ) 


When z -> 0, z °> then 


( II. G. 29) 


SR A ^(E.) 1 / 6 

n \ _ n 1/ 

Sr r=r* 3I/3 r(4/3) 


(II.G.30) 


From (II.G.24) we see that at r = r*, R n vanishes for 
any value of A , therefore 


SR 

W 1 )- -0 


n r=r* 


(II.G.31) 


Differentiating (II.G.28) we have (since ^ = - — ) 

Sr S£ 
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r 2‘Vi""' o/o 

- (-l)" +1 V^ n l/3 J V3 ( -K X nt ) 
Br r-»l L £ 1/2 

_ 2"Vi^ ,/p 

-VF A n 4 J 4/3 (— £ ) 


As 4 approaches zero this becomes 


Br 

7^ } _ 

or r=l 


(-I) n+ \I7* A n 2/3 ( E q ) 1/6 
3 1/3 r(4/3) 


At the outer wall R vanishes but only at A = A 
therefore we must retain both terms of (II.G.15) in 
computing (BR^/Ba^) 


— = 0„ — t C l/2 J,/,] + H„ — [ c l/2 J , /, ] 


+ ^£ U i/ 2j 

3a 1/3 3A 


[ 2 /S ^ 


As C -» 0 the first three terms above vanish but since 


H 2 = 


sin(7 A n - 0 - 5tt/ 12) 
sin( -tt/ 3 ) 


then 


fr* y( - 1 ’ 

V ?/2 
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and 


dR 


n 


)_ 


dA r=l 
n 


(-l) n+1 2^'y 

3 1/b (E„) 1/b r(2/3) 


(II.G.35) 


Now we have 


/ \ (E ) 1/6 3 1/6 r(2/3) 

(c )\ 1} = ^ 

' n 7 i 


(II.G. 36 ) 


7 A 


n 


Similarly 


(C ) (1) = 
N rro 


(_l) n 3 1/6 (E jl/6 r(2/3) 


Vr * 1 Y A 


(II.G.37) 


n 


2 . The fundamental solutions of the second kind 

We must apply the boundary conditions given by (II.F.12) 
Differentiating (II.G.8) we have 


3R, 

Sr 


n 


a. 


Vs J i/3 


2 7 


n 3/2 ) 


- X 2/3 H x J 2/ 3 ( — ^ V) 


sir 


i „3/2' 


(IX.G.38) 


In order for this to vanish as 7 -* 0 (r -* r*), G^ must 
be zero, then from (II.G.12) and (II.G.13) we have 

sin (tt/12 - 0) = 0 


0 = Jo ± nir, n = 0 , 1, 2 


12 


(II.G.39) 


Again n may be taken as zero. We find that 

= 1 

6$ 


H i - 1 


(11. a. 4 o) 



in 


We have a similar condition at the outer wall resulting 


sin(7r/l2 - 7 A n + (f>) = 0 


A n = (n + £•) tt/V, n = 0, 1, 2 ,... (II.G.41) 


This is the equation for the eigenvalues as n becomes 
large. Then, as before. 


H = (-l) n V^ 


(II.G.42) 


For r near r* the solution becomes 


R n^n l/3 ^ ^ 3/2 ) 


-1/3 k 3 


(II.G.43) 


and for r near 1 the solution is 


r = (-i ) n y^ a J 

n ' ' v n 


,3/a, 

-1/3 ( 3 A n ^ ' 


(II. 0.44) 


We will need to get S/Sr (SR n /SA n ) near both walls. 
Differentiating ( II.G.43) we have 


S 

Sr 


SR 


n 


\ SA 


= - A 


1/3 2 

n 3 




n r-»r* 


2\Ie? , 

2z J 0 /i ( 5 z ) 


2Ve7 


+ 2 J ! 2 /3 ( ^ ^ ) 


2/3 k 3 

3 / 2 . 


n 


(II. 0.45) 


As z ■-* 0 (r -» r*), Eq. (II.G. 45 ) vanishes. Near the outer 
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wall we can see in an analogous manner that there Is no 
contribution from the term involving but recall 
that 


~\f r*' sin(ir/l 2 - Y + 0) 
sin( — tt/ 3 ) 


(II.G.46) 


Then 


SG, 

!$k 


(-D n Y 


n 


Vr* j 

VT 


(II.G. 47 ) 


Now 


4 l^s\- (-D n yV? 7 ! 4- [ C 1/2 Ji/s'-r 2 ^ /2)1 

Sr \dA n / ^ Sr /b 0 


2 V/F 1 


As ^ 0 we have 


S ,SR\ (-l ) n+1 YV ^ 7 E 0 l/6 X n 2/3 2 , , , 

I n ' J - 2 n (II. 0.48) 

, 5/0 


Sr \ SA 


n r=l 


3 d/o r ( 4 / 3 ) 


We also need the ordinates at the boundary walls 

,1/3 


R (r*) 


n 


MD = — — I7E 


(e 1 ) 1/S r (2/3) 
(-1 ) n V^ 3 1/3 


n 


(e q ) 1/d r(2/3) 


(n.G.49) 


(11.G.50) 
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Using (II.F.13) and (II.F.14) we may write 


(C ). v 
n' i 


( 2 ) _ 


3 5/6 r(V3) 

2(l-r*)(E i ) I7 ° T X„ 5/3 


(II.G.51) 


n 


(C.) 

' n ' o 


(2) 


(_i) n 3^/6 r(4/3) 

2(1 - r*) E o ) 1//b y 


(II.G.52) 


3. The fundamental solution of the third kind 

f O ) 

For 0^ 1 we have the boundary conditions (II.F.15) 
The condition at the Inner wall yields 


0 = J 75 + mr, n = 0, 1 , 2 , ... 


(II.G.53) 


If we set n equal to 0, we get 


R 


2\Ze7 


n— 


r->r* n 


■ ^- 1/3 2 /2 V 3 <^ ^ 


( II.G.54) 


At the outer wall, by differentiating (II.G.15), we 


have 


sin(j2 - T A n + 0) = 0 


\ = (n + |) Y (11.G.55) 

This Is the equation for the eigenvalues as n becomes 
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large. For r near 1 we have 


2\/ir 


R 


n =(-l)AFC /3 ^ J -l/3 ( ^ X n^ 2 > 

(II.G.56) 


r-»l 


The derivative at the Inner wall Is given by (II.G.30)* 
the ordinate at the outer wall by (II.G.50). We may use 
(II.F.16) to write 


(e ) l/6 r(2/3) 3 1/6 
( C )(3) . . i 


rr i 


(11.G.57) 


y Tv 


n 


For we 

o 


use (II.F.17). At the inner wall 


we 


obtain 


0 = it/ 12 


(II.G.58) 


and 


R n ^ " C /3 zl/2 ^ 


) 


(II.G.59) 


At the outer wall we must have 

sin(y ~K n - 0 - 5ir/l2) = 0 


or 


Tv = (n + 


-) - 

2 Y 


n 


(II.G.60) 


73 



This is the equation of the eigenvalues as n becomes large. 
We now have 


R = (-l) n V^ A 
r-»l 


1/3 4 1/2 


2VE^ 


n 




A n e 3//2 ) (II.G. 61 ) 


The ordinate at the inner wall is given by (II.G.49), 
the derivative at the outer wall by (II.G.33). Then 


(C 

v rr o 


(-i) n 3 1 // 6 (e q ) i/6 r(2/3) 

y 


(II.G.62) 


4. The fundamental solutions of the fourth kind 

For the 15^^ we have the boundary conditions (II.F.19). 
Since (II.F.19) corresponds to (II.F. 17 ) we will have 


A n = (n + (II.G. 63 ) 

as the equation of the eigenvalues with (II.G.49) and (II.G. 33 ) 
giving the inner wall ordinate and the outer wall derivative 
respectively. Then 


(C ) 
n' i 


( 4 ) „ 


,5/6 


’(V3) 

2(1 - r*)(E 1 ) 1/5 y A 


3/3 


(II.G.64) 


For the (II.F. 21 ) must be applied, but this 

correpsonds to (II.F. 15 ). We will therefore have 


*n “ (n + §> ? 


(II.G.65) 
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as n becomes large. The derivative at the Inner wall and 
the ordinate at the outer wall will be given by (II.G. 30 ) 
and (II.G.50). Then 


(C ) 

' n'o 


W = „ 



( II.G.66) 


In practice there will be certain combinations of the 
expansion coefficients and the ordinates or derivatives at 
the boundary wall which contribute to the wall 9's and 
<p l 3 comprising the kernel conditions. 


2(1 - r*)(C n )^ 


( 1 ) 


SR 


Sr 


(!) 2(1 - r*)(E.) l/3 T(2/3) -1/3 

L = A 

1 ' 3 vb r(4/3) 7 

(II.G. 67 ) 


2(1 - r*) (C r ) £ 


m SR (1) 2 (l-r*)(-l) n (E ) 1/6 (E ) l/6 r(2/3) 1 

(1; n ( 1 ) i — A 3 

Sr 


3 1/6 r(4/3) V 


\ 3 


(II.G. 68 ) 


/ x SR ^ 

2(1 - r*)(C (r*) = 

or 


2 (l-r»)(-l)"(B 1 ) l/ 6 (B 0 ) l/6 r(2/3) , , l/3 (ii>0>69) 

~3^T(V3)7V^ n 
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2(1 - r*) (G )(D ^ (1> ( 1 ) = 2(l - r,) ^o )l/3 F(a/3) , -1/3 

n o 8F W F7777T - X n 


3 1/fc> r(4/3) 


.... (II.G. 70 ) 


( 0 n)i 2)R n (2) < r *) - - 


3 1/S r(4/3) _ 5/3 

iTw^TtV 173 r(2/3) y * n 


(C )$ a >R < a )(i) = . J^£.3 l/6 $C r iV 3) 

n 1 n 2( 1 -t>*UE. )^ b T 


(II.G.71) 


•5/3 


2(l-r*)(E 1 ) i/D (E o ) i/D r(2/3) 7 n 


(C )^ 2 ^R ( 2 )(r*) = - 

N rr o n v ' 


(_l) n 3^/ 6 r(4/3) 

2(l-r*) •Vr 5 (E 0 ) 1/6 (E i ) i//6 T(2/3)7 


(II.G.72) 

* -5/3 


n 


(11.G.73) 


( c n ,( a )H n (i, . . _ 37/6 g , ; /3) — , n - 5 /3 (ii . q . 74) 

2(l-r*)(E o ) 1/3 r(2/3) 7 n 


/-> x aR (3 ^ 

2(l-r*)(C )} 3 ' (p*) = 

ar 


2(l-r*)(E 1 ) l/3 T(2/3) _ l/3 

"3 1/6 r(4/3) 7 * n 


(II.G.75) 


<°n/ 3)R n (3) ( 1 ) " * 


(-l) n Vr 7 (E 1 ) l/6 V3~ l _ x 
(E 0 ) l/b 7 An 

76 


(II.G. 76 ) 



(II.G.77) 


(c 

' rr o 


R 


(3) 


n 


(r*) 


(-l) n (E 0 ) l/6 VT _! 

aP y ' X “ 


2 (l-r*)( 0 n )< 3) 


Sr 


(i) 


2 (l-r^)(E o ) l/3 r(2/3 )^ _ l/3 
3 1/b r( 4/3) 7 n 


(c 

4 rr i 


R 


(4) 


n 


(r*) 


(II.G.78) 

3 7/6 r(4/3) ^ _ 5 /3 

2( 1 -r*) (E ± ) 1/3 r(2/3) Y n 

(II.G.79) 


aR ^ 

2(l-r*)(0 )[ 4) — (1) 

or 


(_l) n 3 1 / 2 V ? ? ( Eo ) 1/6 


-1 


(E 1 ) 1/b Y 


n 


(II.G.80) 


m SR 

2 (l-r»)(0 )[ 4) -r 1 (r*) - 

u or 


(-1)" 3 1,/2 ^ _x 

V^OU^Y ^ 


(II. 0.81) 


(o„)i 4 ) R„ ( 4 ) (i) = - 


rr o n 


3 7/6 r(V3) x .5/3 

2(l-r*)(Ej V3 r(2/3) Y " 


(II.G.82) 
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Near the inner wall Eq.(lI.G.7b) is taken as an 
approximation to the actual equation as A n becomes large. 

It is apparent from (lI.G.7a) that if r* is small, A n 
must be very large in order to make (lI.G.7b) a reasonable 
approximation to the actual equation. A comparison of the 
eigenvalues predicted by this method and those obtained from 
the actual solution (see section II. H) shows very good 
agreement for all n > 4 , even if r* is zero. The forms 
(II.G.67) through (II.G. 82 ), however, provide some difficulty, 
particularly those evaluated at the inner wall, or where the 
non-zero boundary condition is applied at the inner wall. 

For r* < 0 . 1 , A^ is not sufficiently large to remove the 
effect of the first derivative term in (H.G.7a). When A n 
becomes very large then (lI.G. 7 a) approaches (lI.G. 7 b) but 
for small r* this value may be so large as to lie outside 
the range of practical interest. What is needed is an 
approximation for the next succeeding constants. Below 
is a logarithmic plot of the quantity 
2(l - r*) (C n ) j 1 ) (dR^ 1 Vdr as a function of A n , 
together with the line representing the approximate 
expression (II.G.67), for r* = 0.02 




o 


* 

u 

1 


OJ 


u 





n 
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Since the constants for n = 1 , 2 , and 3 He on a nearly 
straight line an extrapolation line is constructed to pass 
through the last two. Then the values of the eigenconstants 
are taken as ordinates on this line at the values of 
given by the asymptotic expression. Eventually (i.e. for 
A large enough) the WKBJ expression should become valid, 
and it is used for the region beyond the intersection of 
the extrapolation line and the line corresponding to (II.G. 67 ). 
Both the WKBJ approximation and the extrapolation line have 
ihg form K for all cases. Tables II.G.l through 

II. G. 4 give the values of the constant K for the WKBJ 
approximation and K and S from the extrapolation lines 
for all of the fundamental cases. It can be seen that for 
r* > 0.25 the WKBJ approximation provides an adequate 
expression for the eigenconstants. 

These considerations apply to the absolute values of 
the pertinent functions. The assignment of the correct 
algebraic sign in any case can be made from the expressions 
(II.G. 67 ) to (II.G.82). 

Values of the quantities 7 and tt/y are given in 
Table II. G. 5 and In Figures II.G.l and II. G. 2. 
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TABLE II. 0.1 


ASYMPTOTIC VALUES OF THE FUNCTIONS IN THE FUNDAMENTAL SOLUTIONS OF THE FIRST KIND 



TABLE II. 0.2 


ASYMPTOTIC VALUES OF THE FUNCTIONS IN THE FUNDAMENTAL SOLUTIONS OF THE SECOND KIND 


r* 

(C, 

nW-> 


< C n 

W*> 



'oV r *) 



>oV l > 



K 

K 

s 

K 

K 

S 

K 

K 

S 

K 

K 

S 


WKBJ 

GRAPH 

WKBJ 

GRAPH 

WKBJ 

GRAPH 

WKBJ 

GRAPH 

0 

0 



0 



IS 



2.1389 



0.02 

1.1336 

0.0714 

1.162 

0.22325 

0.0677 

1.462 

11.163 

3.020 

1.423 

2.1389 

2.78 

1.71 

0.05 

1.4039 

0.2065 

1.280 

0.39447 

0.19B 

1.54 

7-8894 

3.86 

1-535 

2.2166 

2.72 

1.70 

0.1 

1.6198 

0.533 

1.442 

0.60112 

0.394 

1.582 

6.0112 

3.94 

1.582 

2.2307 

2.54 

1,6S 

0.2 

1.8261 



0.90404 



4 . 5202 



2.2379 



0.25 

1.6863 

1.122 

1.542 

1 . 0270 

O.863 

1.622 

4.1080 

4.00 

1.66 

2.2360 

2.51 

1.68 

0.3 

1.9335 



1.1384 



3.7947 



2.2341 



0.4 

2.0010 



1.3351 



3-3378 



2.2267 



0.5 

2.047B 

1.78 

1.625 

1.5065 

1.46 

1.645 

3.0130 

3-13 

1.66 

2.2166 

2.45 

1.68 

0.6 

2.0819 



1.6599 



2.7665 



2.2055 



0.7 

2 . 1080 



1.7992 



2.5703 



2.1938 



0.8 

2.1280 



1.9271 



2.4089 



2.1816 



1.0 

2.1483 



2.1403 



2.1483 



2.1463 




TABLE II. 0.3 


ASYMPTOTIC VALUES OF THE FUNCTIONS IN THE FUNDAMENTAL SOLUTIONS OF THE THIRD KIND 


r * 

2(1 - 

r‘)(C n 


(c 

nW 1 ) 


W r *> 


2(1 - 

r-UCnWD 


K 

K 

s 

K 

K S 

K 

K 

S 

K 

K S 


WKBJ 

GRAPH 

WKBJ 

GRAPH 

WKBJ 

GRAPH 

WKBJ 

GRAPH 

0 

« 



0 


00 



4.5477 


0.02 

9 . 9307 

12.40 

0.520 

0.66073 

1.21 1.142 

17.039 

4.38 

0.758 

5.1213 

5-09 0.329 

0.05 

8.3218 

12.12 

0.428 

0.96037 

1.480 1.105 

12.166 

4.72 

0,813 

5-2708 

5-26 0.333 

0.10 

7.4489 

9.84 

0.406 

1.2890 

1.773 1.080 

9 . 3602 

4.91 

0.859 

5.4090 

0.20 

6.8213 



1-7473 


7.1298 



5.5669 


0.25 

6.6623 

7.13 

0.350 

1.9298 

2.190 1.021 

6.5120 

5.33 

0.960 

5.6201 


0.3 

6.5495 



2,0950 


6.0444 



5.6682 


0.4 

6.3898 



2-3855 


5.3599 



5.7422 


0.5 

6.2813 

6.37 

0.338 

2. 6364 

2.86 1.012 

4.8751 

4.60 

0.983 

5.8028 


0.6 

6.2019 



2.8648 


4.5071 



5.8542 


0.7 

6.1398 



3-0703 


4.2154 



5.8998 


0.8 

6.0908 



3.2601 


3-9755 



5-9408 


1.0 

6.0368 



3-6013 


3.6013 



6.0368 



TABLE II. 0.4 


ASYMPTOTIC VALUES OF THE FUNCTIONS IN THE FUNDAMENTAL SOLUTIONS OF THE FOURTH KIND 


r * 

< C n 

W r *> 


2(1 - 

r ->{ c n>l H n 

(X) 

2(1 - 

r*)(c n ; 

l o R A (r ' ) 


W 1 ’ 



K 

K 

s 

K 

K 

s 

K 

K 

S 

K 

K 

S 


WKBJ 

GRAPH 

WKBJ 

GRAPH 

WKBJ 

GRAPH 

WKBJ 

GRAPH 

0 

0 



0 






2.1389 



0.02 

1 . 1336 

0.0709 

1.170 

0.34078 

0.0825 

0.740 

33.037 

65*1 

1.17 

2.1983 

3.70 

1.80 

0.05 

1.4039 

0.2065 

1.282 

0.60830 

0.221 

0.788 

19.207 

27.6 

1.082 

2.2166 

3.27 

1.76 

0.10 

1.6198 

0.417 

1.359 

0.93602 

0.458 

0.836 

12.890 

17.3 

1.07 

2.2307 

2.99 

1.73 

0.2 

1.6261 



1 . 4260 



8.7365 

7.7192 



2.2379 



0.25 

1.8863 

0.976 

1.495 

1.6260 

1.241 

0.939 

9.70 

1.059 

2.2360 

2.63 

1.69 

0.3 

1.9335 



1.8133 



6.9833 



2.2341 



0.4 

2.0010 



2.1440 



5.9638 



2.2267 



0.5 

2.0478 

1.820 

1.635 

2.4376 

2.30 

0.986 

5.2768 

5.80 

1.017 

2.2166 

2.52 

1.68 

0.6 

2.0819 



2.7043 


4.7747 



2.2055 



0.7 

2.1060 



2.9508 



4.3861 



2.1936 



0.8 

2.1280 



3.1804 



4.0751 



2. 1816 



1.0 

2. 1483 



3.6013 



3.6013 



2.1483 
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CONSTANTS IN THE ASYMPTOTIC (LARGE Tv) SOLUTIONS 


r* 

V 

0.0 

0.555360 

0.01 

0.522132 

0.02 

0.516231 

0.03 

0.512245 

0.04 

0.509194 

0.05 

0.506718 

o.q 6 

0.504635 

0.07 

0.502843 

0.08 

0.501275 

0.09 

0.499884 

0.10 

0.498638 

0.20 

0.490731 

0.30 

0.486729 

o. 4 o 

0.484388 

0.50 

0.482937 

0.60 

0.482022 

0.70 

0.481453 

0.80 

0 . 48 lll 8 

1.00 

0.480956 


Table II. G . 5 


7 T/V 

5.65685 
6.OI685 
6.08584 
6.13298 
6.16973 
6.19988 
6 . 22547 
6 . 24766 
6.26720 
6.28464 
6.30034 
6 . 40186 

6.45450 

6.48569 

6.50518 

6.51752 

6.52523 

6.52977 

6.53197 
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CONSTANTS IN THE ASYMPTOTIC CONSTANTS IN THE ASYIVPTOTIC 

(LARGE X) SOLUTIONS (LARGE X) SOLUTIONS 








II. H , Solution of the Sturm-Llouville Equation 


The method used by Graetz and most of the other work- 
ers to generate solutions to the counterpart of equation 
(II. F. 2), for the particular geometry which they considered, 
was the series solution technique of Frobenius. This 
method assumes a solution of the form 


R n< r > * 


00 

z 

m=o 


-m+c 


and by substitution Into the differential equation attempts 
to formulate a recurrence relation among the A n ^ m - The 
principal drawback, for these simpler cases, is the diffi- 
culty with the higher modes of the equation. The A n ^ m 
Involve the eigenvalues A n , and the equation for the 
A Is the condition that the series or its derivative 
vanishes at one of the boundaries. For the higher modes of 

the equation, the A grow very large before diminishing, 

n y m 

and the summing of the series to zero involves the differ- 
ences of large numbers with the result that if accuracy is 
desired many digits must be carried in the computation. 

For equation (II. F. 2) the difficulty is further com- 
pounded by the appearance of the logarithm of the inde- 
pendent variable, r , as a coefficient. To a certain 
extent this difficulty can be circumvented by expanding 
the solution about the point r = 1 , but the resulting 
recurrence relation among the A Is a many term rela- 

II y Hi 

tionship. In fact, the coefficient A_ _ is an explicit 
function of all the preceding A n k = 0,1,2, . . ,m-l, 

which makes practical computation of the solution, at best, 
difficult . 

These difficulties are largely eliminated by the 
iterative method of Berry and de Prima^, This method may 
be essentially summarized as follows: 
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If 


( >h\ 

p 

value, \ , and 


Is 


n 


with 


A 2 - (*?), 


n 'k 

boundary condition 


the kth approximation to the desired 

(R n (r)) k Is a solution to equation (II.P. 

such that (R ), , . . . . 

v n 7 k satisfies the requisite 

only at v = r* ; and if 


2 ) 


1 _ 2 

/ w(r)(R n (r)) k dr = 1 (II.H.1) 


then the next approximation is given by 

(X n> k+ l - <*n>k ± ( R n (1 >>k (H-H.2) 

This sequence of approximations converges monotonically to 
. In (II. H. 2) the plus (+) sign is associated with the 
condition of zero derivative desired at the outer wall 
(r = l) ) and the minus (-) sign with zero ordinate at 
r = 1. 

Computationally, the method consists of assuming a 

value for the slope or the ordinate at r = r* , (whichever 

is not specified as zero by the boundary condition) and 

integrating equation (II.F.2) numerically. The Inner wall 

value Is adjusted to satisfy (II.H.1), then the assumed 
2 

value of is corrected by (II.H.2) and the process is 

repeated. The logarithmic term in the differential equa- 
tion does not pose any fundamental difficulty with this 
method, in fact, w(r) need not be an explicit algebraic 
form at all so that this method can readily be applied to 
these same problems In turbulent flow. 

The computation of these solutions has been performed 
on the Burroughs 220 Electronic Digital Computer at the 
Stanford computation center. The nature of the computation 
is such that the values of the expansion coefficients are 
obtained directly from (II.F.7) using a simple numerical 
integration. The condition ( II.H.1) makes the denominator 
of (II.F.7) always equal to unity. Details of the compu- 
tational procedure are contained in appendix E. 
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The eigenvalues and the pertinent combinations of con- 
stants are given in tables II.H.l to II. H. 4. The same 
quantities are presented in graphical form in figures 
II.H.l to II.H.20 to facilitate interpolation for the inter- 
mediate values of radius ratio. 

The values for parallel planes (r* = l) were provided 

by McCuen 11 , those for the circular tube (r* = 0) are from 
> 20 
Lipkis 10 and Siegel, Sparrow, and Hallman 

For the two limiting geometries, the tube and the 
parallel plane channel, the constant N (which is a func- 
tion of r* ) has the values 1/2 and 3/32 respectively if 
x and r are defined in the manner used here. It seems 
customary* to absorb that constant Into the eigenvalue 
resulting In values which are somewhat different than those 
reported here. The practice would produce some difficulty 
for the case at hand since N — > °° as r* — > 1 , which 
would result in A becoming very large. To preserve the 
continuity of the eigenvalues as r* — ^ 1 this constant 
Is excluded from A . To transform the previously reported 
eigenvalues for the tube and the parallel plane channel 
into the values reported here requires multiplication by 
y/2 and \/32/3 respectively. 

An analogous difficulty arises In regard to those 
functions which Involve the derivatives at the walls. 

SR n 

C — as r* — > 1 

n a— 
or 


but the product 


2(1 - r*) C n 



* IQ 

C.f. Sellars, et . al. 
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remains finite. Further, it is this second form which is 
desired in computing the $ T s since this is equivalent to 

3R n 

D h C n 3F 1 

All the functions involving the derivatives are reported 
in the above form. 
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TABLE II.H.l 


FUNCTIONS IN THE FUNDAMENTAL SOLUTIONS OF THE FIRST KIND 




2(1 - r*) 

<v (1) 

<*n> (1) 

•< C n>i 1,R A< r ' 

4.74805 

10.8763 

16.9948 

23.1043 

22.5440 

118.294 

228.824 

533.610 

-10.2263 

- 6.39215 

- 5.06568 
' 4.35330 

4.93941 

11.1805 

17.4093 

23.6290 

24.3978 

125.003 

303.085 

538.332 

- 6.69710 

- 4.45700 

- 3.64722 

- 3.20042 

5.10557 

11.4367 

17-7580 

24.0720 

26.0669 

130.799 

315.346 

579.459 

- 5.17154 

- 3.61850 

- 3.02917 

- 2.69415 

5.32304 
11.7580 
18. 1952 
24.6305 

28.3348 

138.252 

331.065 

606 . 663 

- 4.05663 

- 3-00828 

- 2.51348 

- 2.31513 

5.44515 

11.9303 

18.4308 

24.9339 

29.6562 

142.332 

339.693 

621.698 

- 3.63663 

- 2.77072 

- 2.38880 
- 2.15661 

5.49209 
11.9937 
18. 5177 
25.0463 

30. 1630 
143.848 
342.906 
627.316 

- 3.43234 

- 2.60888 

- 2.27786 

- 2.05809 


2(1 - r* ) 

•< c n>i 1, "A< 1 > 

+ 0.80760 

- 0.549498 

+0.452437 

-0.397809 

+1.03835 

-0.730683 

+0.611777 

-0.543755 

+1.29745 

-0.938758 

+0.795670 

-0.712271 

+1.82903 

-1.37025 

+1.17600 

-1.05935 

+2.46559 

-1.88131 

+1.62265 

-1.46516 

+3.43234 

- 2.60888 

+2.27786 

-2.05809 


•(Cn)' 1 ^*) 

-40.2547 

+27.3773 

- 22.5521 

+19.8221 

-20.7790 
+14.6215 
-12.2431 
+10.8819 
-12,9769 
+ 9.38927 

- 7.95812 
+ 7.12413 

- 7.31778 
+ 5-48090 

- 4.70388 
+ 4.23720 

- 4.93127 
+ 3.76263 

- 3.24530 
+ 2.93030 

- 3.43234 

+ 2.60888 

- 2.27786 

+ 2.05809 


+3.17904 

+2.35347 

+2.01421 

+1.81137 

+3.22168 

+2.39706 

+2.05363 

+1.84885 

+3.25568 

+2.43589 

+2.09036 

+1.88347 

+3.29990 

+2.49651 

+2.14955 

+1.93866 

+3.34333 

+2.55482 

+2.20445 

+1.99079 

+3.43234 

+2.60888 

+2.27786 

+2.05809 


FUNCTIONS IN THE FUNDAMENTAL SOLUTIONS OF THE SECOND KIND 


r* 

n 

(X )< 2 > 
' rv 

<# (2> 

(C n )< 2) R n (r*) 

0 

1 

7 . 16654 

51.3593 

0 


2 

12,9508 

167.724 

0 


3 

I8.6636 

348.329 

0 


4 

24.3532 

593. 080 

0 

0.01 

1 

7.50711 

56.3567 

-0,00347335 


2 

13.6217 

185.551 

-0.00178909 

0.02 

1 

7.53889 

56.8349 

-0.00665350 


2 

13.7049 

187-825 

-0.0033414 


3 

19.807B 

393-348 

-0.0021627 


4 

25.8930 

670.448 

-0.001603 

0.035 

1 

7.55723 

57.1117 

-O.OIIO6II 

0.05 

1 

7.56239 

57.1897 

-0.0151141 

2 

13.8132 

190.803 

-0.0071613 


3 

20.0134 

400.537 

-0. 0044722 


4 

26 . 2024 

686 . 569 

-0.0031492 

0.075 

1 

7.55942 

57.1448 

-0. 0212168 

0.1 

1 

7.55041 

57.0087 

-0.0266349 


2 

13.8836 

192.753 

-0.0117253 


3 

20.1791 

407.198 

-0.0069693 


4 

26.4683 

700.572 

-0.0047401 


5 

32.7568 

1073.01 

-0.0034799 


6 

39.0459 

1524.58 

-0.0026878 

0.175 

1 

7.51572 

56.4860 

-0.0397586 

0.25 

1 

7.48406 

56.0112 

-0.049403 

2 

13.9539 

194.710 

-0.019141 


3 

20.3940 

415.917 

-0,010610 


4 

26 . 8276 

719.722 

-0.006890 

0.375 

1 

7.44514 

55.4301 

-0.060955 

0.5 

1 

7.42034 

55.0615 

-0.068661 

2 

? 

13.9888 

20.5195 

27.0385 

195.686 

421.048 

731.082 

-0.024215 

-0.012888 

-0.008178 

0.7 

1 

7.39945 

54.7518 

-0.076527 

1.0 

1 

7.39139 

54,6329 

-0.003014 

2 

14.0023 

196.064 

-0.027778 


3 

20.5686 

423.067 

-0.014468 


4 

27 . 1205 

735.519 

-0.008670 


(c n )< 2) VU 


+0 . 00179521 
-0.00075895 

+0. 0035018 
-0.0014612 
+0.0008562 
-0.0005839 
+0.0059587 
+0.0083162 
-0.0033617 
+0.0019246 
-0.0012802 

+0.0120597 

+0.0156020 

-0.0060579 

+0.0033763 

- 0.0022056 

+0.0015744 

-0.0011910 

+0.0252307 

+0.033643 

-0.012135 

+ 0.006510 

-0.004146 

-HD. 045580 

+0.055527 

-0.018997 

+0.009971 

- 0.006276 

+0.068417 

+0.083014 

-0,027778 

+0.014468 

- 0.008670 


( c n>rv-*> 

+0.201742 

-0.087555 

+0,052797 

-0.03664 

+0.179521 

-0.075893 

+0,175092 

-0.073059 

+0.042809 

-0.029010 

+0. 170252 

+0.166324 

-0.067237 

+0.038488 

-0.025607 

+0.i60796 

+ 0.156018 
- 0.060580 
-HD. 033790 
- 0.022058 
+0.015740 
-0.011912 

+0.144177 

+0.134564 

-0.048534 

+0.026042 

-0.016590 

+0.121550 

+0.111053 

-0.037994 

+0.019945 

-0.012553 

+0.097750 

+0.083014 

-0.027778 

+0.014468 

- 0.008670 
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TABLE II. H. 3 




FUNCTIONS IN THE FUNDAMENTAL SOLUTIONS 

OF THE THIRD 

AND FOURTH KINDS 






2(1 - r*) 


2(1 - r*) 


r* 

n 


>£ 4! <*S> 

■< C n>i 3)R A< r *> 

(C n )'3) Rn(l) 

•< C n>o 4)R A< r# > 

(C )i 4 ^R (l) 
' rro n' ' 

0.02 

0 

1.52956 

2.53626 

-31.7792 

-1.07360 

-53-5098 

-1.80732 


1 

8.45591 

71.5024 

- 7.59043 

+0.108022 

+ 5.36368 

-0.076332 


2 

14.7044 

216.219 

- 5.52771 

-0,056343 

- 2,81542 

-0.028697 


3 

20.8730 

435.682 

- 4.60870 

+0,037862 

+ 1.87650 

-0,015416 

0.05 

0 

1.82368 

3. 32580 

-17.1123 

-1.09236 

-21.8599 

-1.39542 


1 

8.73469 

76,2948 

- 5.17349 

+0.136946 

+ 2,73988 

-0.072526 


2 

15.0887 

227 , 668 

- 3.93218 

-0.073469 

- 1.47056 

-0.027476 


3 

21.3656 

456.490 

- 3.36102 

+0,050123 

+ 1.00256 

-0.014951 

0.1 

0 

2.05027 

4.20362 

-11.2121 

-1.11276 

-11.1298 

-1.10459 


l 

8.97827 

80.6094 

- 4.12691 

40.168675 

+ 1.68692 

-0.068948 


2 

15.4168 

237.678 

- 3.23924 

-0.092579 

- 0.924952 

-0.026464 


3 

21.7847 

474.573 

- 2.81551 

40.063855 

+ 0.638696 

-0.014486 

0.25 

0 

2.42830 

5*89665 

- 6.98949 

-1.15233 

- 4.60952 

-0.759954 


1 

9.31079 

86.6908 

- 3.35613 

40.230576 

+ 0.922335 

-O.O63367 


2 

15.8447 

251.055 

- 2.72768 

-0.130316 

- 0.521331 

-0.024907 


3 

22.3254 

498.423 

- 2.40759 

40.090973 

+ 0.363942 

-0.013752 

0.5 

0 

2.76602 

7. 65088 

- 5.30242 

-1.19491 

- 2.38986 

-0.538562 


1 

9.53544 

90.9247 

- 3.04516 

+0.297663 

+ 0.595329 

-0.058182 


2 

16 . 1056 

259.389 

- 2.51847 

-O.171774 

- 0.343531 

-0.023431 


3 

22.6421 

512.665 

- 2.23664 

+0.120832 

+ 0.241646 

-0.013055 

1.0 

0 

3.11796 

9.72165 

- 4.35312 

-1.24042 

- 1.248436 

-0.358038 


l 

9.71413 

94.3647 

- 2.854332 

+0.383198 

+ 0.393268 

-0.051454 


2 

16.2682 

264.654 

- 2.387368 

-0.226323 

- 0.226400 

-0.021463 


3 

22.8118 

520.379 

- 2.127292 

+0.160551 

+ 0.160704 

-0.012129 






TABLE II. H. 
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FUNCTIONS 

IN THE FUNDAMENTAL SOLUTIONS 

OF THE THIRD 

AND FOURTH KINDS 







2(c - r*) 


2(1 - r*) 

r* 

n 


< C n>i 4 ) V r *> 

•< c n >i 4 )R A(i> 

< C n>o 3 ) V r *> 

•<°n>P )R A< 1 > 

0 

0 

3.824549 

14.62718 

0 

0 

-1.46622 

+2.99516 


1 

9.445581 

89.21900 

0 

0 

+0.80248 

+2.17696 


2 

15.09444 

227.8421 

0 

0 

-0.58709 

+1.85152 


3 

20.74805 

430.4816 

0 

0 

+0.47490 

+1.66072 

0.02 

0 

3.95731 

15.6603 

-0.0127464 

40.0289949 

-1.44975 

43.29778 


1 

9.94704 

98.9435 

-0.0047162 

-0.0150085 

+0.750414 

42.38808 


2 

15.9795 

255-346 

-0.0027714 

+0.0106132 

-0,530669 

42.03225 


3 

22.0248 

485.094 

-0.0019160 

-0,0083610 

+0.418050 

41.82426 

0.05 

0 

3.93106 

15.4532 

-0.0305157 

40.0716530 

-1.43305 

43.36492 


1 

9.99141 

99.8283 

-0.0104624 

-0,0357396 

+0.714786 

42,44169 


2 

16.1156 

259.711 

-0.0058357 

40.0246326 

-0.492656 

42.07949 


3 

22.2616 

495.518 

-0.0038742 

-0.0190190 

+0.380383 

41,86734 

0.1 

0 

3.86760 

14.9583 

-0.057923 

+0.141044 

-1.41047 

43.43456 


l 

9.99430 

99.8861 

-0.017888 

-0.066845 

+0.668553 

42.49829 


2 

16.2096 

262.750 

-0,009377 

+0.044671 

-0.446729 

42.12819 


3 

22.4544 

504,200 

-0.005966 

-0.033772 

+0.337525 

+1,91074 

0.25 

0 

3.68006 

13.5428 

-0.128700 

+0.340089 

-1.36035 

+3.59471 


1 

9.93693 

'98.7426 

-0.031544 

-0.143358 

40.573406 

42.60594 


2 

16.3028 

265.783 

-0.014912 

40.090833 

-0.363314 

42.21309 


3 

22.6973 

515.170 

-0.008924 

-0.066513 

40,266097 

41.98324 

0.5 

0 

3.43679 

11.8115 

-0.221934 

40.653514 

-1.30702 

43.84871 


1 

9.84575 

96.9388 

-0.042427 

-0.240047 

+0.480096 

42.71634 


2 

16.3249 

266 . 502 

-0.018654 

40.146212 

-0.292413 

42.29196 


3 

22.8169 

520.611 

-0.010781 

-0.105092 

4<>. 210151 

42.04857 

1.0 

0 

3.11796 

9.72165 

-0.358038 

+1.248436 

-1.24842 

44.35312 


1 

9.71413 

94.3647 

-0.051454 

-0.393268 

+0.383198 

42.854322 


2 

16 . 2682 

264.654 

-0.021463 

+0.226400 

-0.226323 

42.387368 


3 

22.8118 

520.379 

-0.012129 

-0.160704 

+0. 160551 

+2. 127292 
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II . I . The Fundamental Solutions 


We are now prepared to construct the various funda- 
mental solutions using the relationship 


0 (k) = ^(k) 0 (>)) (II. I. 1) 

j 3 J fd 

and the series solution for 0 ^) , (II.F. 8 ). The fully 
developed solutions are taken from section II. D. 

The pertinent non-dimensional quantities (as outlined 
In section I.B) may be readily obtained from (II. I. l) and 
the defining equations. Let us formulate some general 
relationships of interest before turning to the specific 
cases 


<4^ (x) = -2(1- r* ) 


( x ) = 2 (1 - r* ) 


80 


00 


8 r I—’ „ 

r=r* 


80 


(k) 


i 8r 


(II. I. 2) 


(II. 1.3) 


r=l 


The 0 ^) may be obtained from (I.B. 11) or, since 
/. \ m J _ 

= 0 at x = 0 , from consideration of an energy 
mj 

balance on the flowing fluid. We see that 


m J - = ( r* 4^ (x) + 4^ (x)) (II. 1.4) 

87 1 + r* i J OJ 


or 

0^(7) = ( ) / (r*4^(x) + dx (II. 1.5) 

m J 1 + r * o J J 

The conventional treatment of problems of this sort gener- 
ally involves the Nusselt Modulus. It Is related to the 
0 's and the 3 >'s by 
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00 

Nu „ j 


1. The first kind 


AV - 9 J£'] 


(II. I. 6) 


-<V 8 > (1 * 


+ 2 {C n )[^ R^(r) e n 


= ♦I?')) 


ii 'fd 


°° ( , v /dR'^M 

- 2 2(1 - r * ) ( C ) j 1 —2— 

^ L1 -L jTT . “ 


-(*n) (1) * 


dr / r=r* 


oi ^ ' ol 'fd 


+ 2 2(1 - r*)(C ) 

n=o 


m /dR^h 

) (1) - n e 

n 'i ,— — 

I dr ' r=l 


-K 2 > (1) * 


0^)( x ) _ g(l)\ + 

mi ' ' mi /fd + 7T 


7¥T 2 2(1 - r*)(C„) (l) 


dr L =r , 1 dr L J 


-(^) (l) x 


(1) 

Nu ii (x) - ; — xttt=; 


1 - 


(II.I.ll) 


100 



(II. I. 12) 


m - 

Ni/ } ( x ) = — t i v _ 
01 


= 1.0 


e^^x) = 0.0 


0o 1 ) (r,x) = 1 - fr^ + S 

u n=o 


(II. I 
(II. I 

-(^) (1) K 

(II. I 


^ ) (-) = ^o ) )fd + S 2(! - r *)( c n >o 


( 1 ) 


00 


n=o 


/ dR 


(1) 


n 


dr 


-(^) (1) X 


(II.] 


r=l 


•ii’t*) - - j 0 2 < 1 - 


I dR 


( 1 ) 


n 


x 


(11.: 


dr 


r=r* 


0 ( 1 )(7) = + jy-J— 5 -y 2 2(1 - r*)(C ) K 

mo ' ' mo 'fd (1 + r* ) n _ Q rJ u 


( 1 ) 


/dR. 


( 1 ) 


r* 


n 


fdR, 


( 1 ) 


n 


-(A 2 ) (l) 7 
v n * 




dr / r=r * \ dr I r^l J 

(11. 


.13) 

• 14) 

• 15) 

:.i6) 

: - 17 ) 


1.18) 
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(II. I. 19) 


NuO>m. 

00 1 - 8 ^(x) 

mo x ' 

Nu^)(x) . - 4°| ( ^ 
mo x 7 

0 (l) (x) = 1.0 

00 v ' 

e^^x) = o.o 


2. The second kind 


e[ 2 )(r,x) = — 2 

1 / "1 « -H. T Id. \ 


(1 + r*)M(l - r* )' 


r 4 , Br 2 - ( r* ) 2 M , ( 1 + in r* 

T + x“ r ' r + f — 


00 (?) (?) -(A 2 ) (2 >x 

+ 2 (C n )[ 2 ) R^ 2 )(7) e n ( 

n=l 1 n 


(2) _ 0 (2) n . 

11 ml '' fd + 1 + r* X 


+ 2 
n=l 


( c n )i 2) 4 2) (^*) e 


-(X n 2 ) (2 >x 


( S (2) - e(2)) + in 

' oi ml ; fd 1 + r* 


+ 2 (C n )* 2) R* 2) (l) e 


-(x 2 ) ( % 


, 4r* \ — 

(t~ o. t.« ) x 
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( II . I • 27 ) 


$( 2 ) (x ) = 1.0 

$ol^(x) = 0.0 


Nu 


( 2 ) 

ii 


(x) 



(11. 1. 28) 

(11. 1. 29) 




(11. 1. 30) 

(11. 1. 31) 

(11. 1. 32 ) 

(11. 1. 33) 

(11.1.34) 

(11. 1. 35) 
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(11.1.36) 


jj u ( 2 ) _ 1 ♦ 0 

00 e^ ; (x) - 


mo 


3 . The third kind 


ep>( 7 ,x)= 


1+2 

n=o 


< C A 


< 3 ) (R n (r))j 3) 


-(* 2 )| 3 >x 

' n 'i 


(II. I. 37) 


♦ii^x) = - 2 2(1 - r*)(C ) ( 

V'i /-v 1 -* _L 


o) M ( v r »i 3) r 


n=o 


dr 


r=r* 


-^>i 3) * 


= 0.0 


®o 3 ) (x) 

e< 3) (?) = 1 + Z ( c n )i 3 '0Ul)), l3J e 

n=o 1 


( 3 )„ nu( 3 ) - (X n)i 3) * 


si 3) (*) 

e i? 4 x) 


= 1.0 


( 11 . 1 . 38 ) 

(IX. 1.39) 

(11. 1.40) 

(11. 1.41) 


= l 


4 r* 


+ TT-pr^ - r‘)(Cn)} 


(3) 


k(r))p> 




dr 


r=I '* Un)P 


Nu^^(x) 
11 v ' 


4 >[^(x) 

1 - 


( 11 . 1 . 42 ) 


( 11 . 1 . 43 ) 


For this case the limiting value of the Nusselt number as 
x becomes large Is of considerable Interest and was not 
apparent from the asymptotic analysis of section II. D. 

If we write out the series expressions appearing in 
(II.I. 41 ), and notice that as x -> oo only the leading 
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terms of these series make any contribution, we obtain 


Nu (3)\ _ 1 + r * (a 2 ) (3) 

N ii ' fd ^Fr* ' A l'i 


(II. I. 44) 


4 3 > ( 7 ,x) = 1+2 ( 0 n )< > 3 ) (R n (r))< 3) 


v n 'o 


n=o 


e' 3 ) (x) = l.i 


f4 3) (x) = 1+2 (C n )p ) (R n (r*))J, 3) e 


,(3) --<^ 3fc 


. . . (ii. I. 45) 

(11. 1. 46) 

(11.1.47) 


n=o 


n=o 


n 'o 


*< 3) (x) = 2 2(1 - r*)(C n ) 

n=o 

-(A?) 1 3 >x 


(3) K( R n( r )>o 3) ] 


dr 


( r=l 


n ' o 


4o } ( x ) = 0,0 


(11. 1. 48) 

(11.1.49) 


0 ^ 3 ^(x) = 1 - 
mo v 1 


ITT-FTT Jo 2(1 ' r * )( 0 n ) c 3> 


fd[(R n (r))< 3) ] \ 


-< X n)o 3>x 


dr 


/-i IS W 


, v 3 ^ ( x ) 

Nu^ 3 ^(x") = °°Mj 

oo v 1 - e^ 33 (x) 

. mo ' ' 


(II. I. 50) 


(II. I. 51) 


(II. I. 52) 


The solution of the third kind, when the temperature step 
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is applied at the outer wall is the extension of the clas- 
sical Graetz problem and reduces to it at r* = 0. 

4. The fourth kind 




ii 


- St"- r*) *" 7 + n f 0 (°n)[ 4) ( R n( 7 ))| 4) 

e -(^ 

(II. 1.53) 


-<^>i 4) * 


. . . . (11.1.54) 

0.0 

(II. I. 55) 

1.0 

(11.1.56) 

^oi ) 4d + J 0 2(1 - r *>< C n>i 4) 


( d [( R n ( 7 ))l 4> ]\ c - (X n)i 4)5r 

(II. I. 57) 


dr 


- e iV) 


_ i 2 2(1 - r* ) (C 

mi 1 + r* ' ''^n'i 

n=o 



.-<*>i 4 * 


NU “ > " ^ 


Nu 


(4) _ *ol ) ( x ) 
01 “ 


(11.1.58) 


(II. I. 59) 


(II.I.60) 
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e^tr.x) • 2(1 - r* T ln + n ^ 0 ( C r7o ^ R n^ r ^o 


-U 2 ) (4) x 

v n 'o 


(n 


0 ( 4 ) ( 7 ) = 
00 v ' 


+ j o ( 0 n ) i 4 )(R n ( 1 ) )i 4)e " Un) ° X (I1 


- 


0.0 


(II 


oo ' 7 


1.0 


(II 


^ 4 ) (x) = 


^ 4) ) fd - S 2(1 - r *)(0 n )( 4) 

ia n=o 


[( ( r ) ) o 






dr 


r=r* 


(II 


= * 2 d--)(c n )i 


(4)x . 4r* Z o/n „*\f n ^*0 


n=o 


(RntO)! 4 ’] 


-t"n)o 4)? 


Nu'^lx) = 


Nu^t*) = 


dr 

1.0 


r=r * (*n } o 


^<*> ‘ 0 mo j ( x ^ 


(II 

(II 

(II 


mo 


1.61) 

1 . 62 ) 

1.63) 

1.64) 

.1.65) 

. 1 . 66 ) 

• 1 . 67 ) 

. 1 . 68 ) 
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The values calculated for the various geometries are 
presented in tables II. I. 1 to II. I. 8 and figures II. I. 1 
to II. I. 12. 

As is indicated by the preceding formulas, the values 
of interest at any axial station, x , are obtained by 
adding or subtracting the series of the eigenfunctions to 


the fully established or terminal values. In all the cases 
where x 1.0 the terminal values have (or very nearly 
have) been obtained. As x diminishes the series repre- 
sents a progressively greater part of the local value. 

Since the series contains an exponential term. 


e , the larger values of x receive only the contri- 

bution of the first few terms, which are known very pre- 
cisely from the machine calculation. For values of 
x = ^ 0.01 , however, the higher terms (those beyond 
n = 3 or 4 in the present case) have begun to contribute. 
These terms are obtained from the graphical method 


described in section II, G and from the asymptotic (large A) 
solutions. While the machine calculation easily provides 
5 or 6 significant figures, the terms beyond the third or 
fourth are known only to 1 or 2 percent and, accordingly, 
the level of confidence in the values of the fundamental 
solutions diminishes with x . 

Those quantities which start from zero at 5c = 0 and 
Increase to their terminal values somewhere downstream are 
particularly affected by this uncertainty since, for small 
values of x , the absolute values of these quantities are 
provided by the differences of nearly equal quantities and 
the uncertainty in the series represents a- large fraction 
of the residual values. 

An effort is made, in the tables, to represent the 
confidence In the recorded values by the number of signi- 
ficant figures retained. Fortunately the physical nature 
of the problem provides some internal checks on the 
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accuracy, i.e. for small x approaches , 9 ^ 

approaches 0^^ * e ^ c> From these considerations it 
appears that the values of the series used for the wall 
values at x = 10 - ^ are within + 1 percent of the actual 
or correct values. At this value of x , ^ and 

( p) ( 4 ) 

are reliable to ± 1 percent, and 0^' and 0^ are 
(because they have diminished in absolute value) reliable 


to + 5 percent . 

The opposite wall values (e.g. 0 ^ ) diminish rather 

rapidly and in all the cases for x < 0.01 they are sub- 
stantially zero (less than 0.1 percent of terminal value). 

For very small values of x the Leveque approxima- 
tions provide a good estimate for the applicable quantities. 
For all cases where the non-zero boundary condition is at 
the outer wall, and for the inner wall heated cases with 
r* >. 0.5 , the Leveque solution agrees well with the 
computation at x = 10 ^ . For small values of r* , when 
the heating is at the inner wall it appears that the 
approximate solution will not become valid until x is 
very small. For r* = 0.02, x 10 seems to be required, 
but for r* > 0.05, x 1 10”^ appears adequate. 

This difficulty arises from the neglect of the curva- 
ture term in equation (IX.E.5). For small values of r* 
the depth of penetration of tne temperature signal must be 
very small before this is justified. 

The graphical representation of the fundamental solu- 
tions provides a very convenient method of superposition 
to construct the more complicated boundary conditions. 
Examples of this are presented in appendix G. 
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THE FUNDAMENTAL SOLUTIONS OF THE FIRST KIND 
TEMPERATURE STEP AT THE INNER WALL 
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FUNDAMENTAL SOLUTIONS OF THE FIRST KIND 
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THE FUNDAMENTAL SOLUTIONS OF THE SECOND KIND 
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THE FUNDAMENTAL SOLUTIONS OF THE THIRD KIND 
TEMPERATURE STEP AT THE OUTER WALL 
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FUNDAMENTAL SOLUTIONS OF THE FOURTH KIND 
INNER WALL HEATED 
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the fundamental solutions of the fourth kind 

OUTER WALL HEATED 
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THE FUNDAMENTAL SOLUTIONS OF THE SECOND KIND 
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THE FUNDAMENTAL SOLUTIONS OF THE SECOND KIND 
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THE FUNDAMENTAL SOLUTIONS OF THE THIRD KIND 
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FIGURE 0.1.7 


THE FUNDAMENTAL SOLUTIONS OF THE THIRD KIND 
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FIGURE n.1.8 
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III. EXPERIMENTAL WORK 


III. A. Description of Apparatus . 

The apparatus for the present investigation consists of 
a positive pressure system using air as the working fluid. 
The boundary walls of the annular passage are inconel tubing 
which themselves form the resistance heating elements. 
Electrical power is supplied by a two channel regulated AC 
power supply, each channel having 1 KVA capacity at 200 
amperes maximum current. Figure III. A. 2 is a diagram of 
the power supply circuit including the principal instrumen- 
tation . 

The test section is 6 feet long, although in practice 
only 4 feet of this is heated. In this series of tests the 
initial 2 feet was used as a "starting" length to allow the 
velocity profile to develop. Two sizes of outer tube are 
used with inside diameters of approximately 1.96 and 1.00 
inches. These, together with core tubes with outside 
diameters of approximately 0.375 and 0.500 inches, provide 
four radius ratios of 0.191, O. 256 , O. 376 , and 0-501. 

Considerable care was taken to secure close tolerance 
injconel tubing. The specifications on the tubing are: 

a. 2 inch outer tube ; 

I.D., 1.96 inches ± 0.0025. 

wall thickness, 0.020 inches ± 0 . 001 , 
ovality, ± 0.010 inches, 
straightness, ± 0.010 inches/foot; 

b. 1 inch outer tube, 0.375 and 0.500 inch core tubes 
I.D. or O.D., specified values ± 0.001 inches, 

wall thickness, 0.020 inches ± 0 . 001 , 
ovality, ± 0.002 inches, 
straightness ± 0.002 inches/foot. 

Inconel was selected because of its high specific resistivity 
( 98 - 1 - cm) which enables the use of relatively low heating 

currents . 
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In view of the close tolerance on the wall thickness 
(the pieces actually used have variations of wall thickness 
of less than ± 0.0003 Inches) electrical heating provides 
essentially constant heat flux at the walls . Accordingly , 
the experimental apparatus provides values for the funda- 
mental solutions of the second kind (and superpositions 
thereof) directly. 

The assembled test section Is mounted, vertically, on 
an 8 foot piece of aluminum channel 6 inches wide. The 
outer tube is supported In three lucite holders, two at the 
extreme ends of the tube and one at the point of the electrl 
cal connection to the heated tube. The core tube is held by 
two aluminum holders mounted 8 inches beyond the ends of the 
outer tube. To minimize eccentricity the holders were 
clamped together and bored simultaneously. While clamped, 
two aligning holes were drilled and reamed through all of 
the pieces. During the assembly on the channel these 
aligning holes were brought into line with a transit equipped 
with a traversing rack. Machined brass plugs with a central 
hole made with a #67 (0.032 inches) drill were pressed into 
the aligning holes. These holes were readily distinguished 
by the transit during the assembly, accordingly, it is felt 
that the axial misalignment of the holders Is less than 
± 0.015 inches . 

A positioning bracket that contains the center tube fits 
Into a bored recess in the end holders and bears against an 
" 0 " ring. A teflon plug fits over the tube and bears against 
it and the Inside of a tapered hole in the positioning bracket 
effecting a seal and centering the tube with respect to the 
bracket. When installed, the core tube is put in tension to 
straighten any residual bow. 

Four brass aligning spheres were made each having a 
diameter 0.010 inches less than the annular gap provided by 
the four geometries. During the running of the tests the 
concentricity is periodically checked by lowering the 
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appropriate sphere into the test section from the top. This 
insures that the axial misalignment will be less than 4$ of 
the smallest gap width under all conditions. 

The air enters at the lower end of the channel which 
is sealed to provide a plenum chamber 6 inches on a side. 

To facilitate the right angle turn made by the flow 
a set of sheet metal turning vanes is mounted inside the 
plenum chamber. The entrance is a wooden nozzle with the 
converging wall profile a quadrant of an ellipse. The up- 
stream face of the nozzle seals against an "0" ring in the 
top wall of the plenum. The joint between the nozzle and the 
tube section was carefully filled, sanded, and varnished to 
provide a continuous wall surface and to avoid tripping the 
boundary layer during the laminar flow runs. Details of the 
assembly are shown in Fig. III. A. 3 . 

Just upstream of the nozzle, three screens are placed 
to reduce the turbulence level entering the test section and 
promote even distribution of the air flow. With both the 
inner and the outer tubes removed, a hot film anemometer 
(Ref. 17) was used to measure the turbulence intensity at the 
center of the exit plane of the nozzle having the 1.96 inch 
throat. The results are as follows: 


Turbulence 

Condition Velocity-f t/sec Intensity- u 1 /U,$ 


No screens 

c 

50 

— L 11 

6.9 

One fine screen 

23 

1.3 


64 

3-2 

Three screens 

10 

0.30 


46 

0.37 


79 

O .58 
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The velocity profile at the nozzle exit plane was surveyed with 
a Pitot-static probe to insure that no serious maldistribution 
of flow resulted from the rather short radius turning of the air 
stream. The results of a survey are shown in Fig. III.A.l. 
These profiles were measured along a diameter of the nozzle 
exit normal to the mounting channel (i.e. along a radius of 
the curved streamlines) as the skewness of the flow is 
greatest along that line. [ 



FIG. III.A.l 

The air is supplied to the test section by a centrifugal 
blower producing a maximum static head of 30 inches of water. 
The air supply is shared with the companion apparatus, the 
parallel plane channel as described by McCuen . The air 
is metered by flat plate orifices. Either of two parallel 
metering sections is used, the sizes are as follows: 


Duct diameter 

Orifice diameter 

P = d/ D - 

6.367 

4.225 

0.664 

6.367 

1.910 

0.300 

1.498 

0.824 

0.550 

1.498 

0.375 

0.250 


v 
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A diagram of the metering system is shown in Fig. III. A. 4. 
System static pressure is measured by a 30 inch vertical 
water manometer, differential pressure by one of three 
parallel manometers, a l6 inch vertical water manometer, or 
inclined oil manometers having ranges of 1 and 5 Inches of 
water. The differential manometers were adjusted by 
comparison with a micro-manometer and the entire metering 
system was run in series with portable laboratory nozzle 
tanks whose calibration is regarded as well established. 

The total range of air flow which can be measured is from 
1 to 750 cfm with a probable uncertainty of about ± 1 $. 

The power generated in the test section is determined 
by an ammeter in series with and a voltmeter across the re- 
sistive load. The ammeter used is a Weston model 901, the 
voltmeter a Hewlett-Packard model 400H, VTVM. A switching 
arrangement is used to permit the use of only one ammeter 
and one voltmeter even when both tubes are heated (see 
Fig. III. A. 2). 

The fluid mean temperature is determined as a function 
of axial position by adding the temperature rise due to the 
energy input to the inlet temperature. The inlet temperature 
Is measured by thermocouples fastened to the turning vanes 
just upstream of the nozzle entrance. 

At the outer wall, temperatures are measured by lron- 
constantan thermocouples spot welded directly to the outer 
surface of the tube. When the outer tube is heated a correc- 
tion is made for the deviation of the outer surface tempera- 
ture from that of the Inner (see appendix A.). Power Is 
supplied to the outer tube through OFHC copper rings Soft- 
soldered to the outside surface. Thermocouple locations 
together with the positions of the current rings are shown 
in Fig. III. A. 5. 

In the core tube, thermocouples are placed as shown 
in Fig. III. A. 6. When Installed, the core tubes are position- 
ed so that the thermocouples approximately oppose those on 

* 
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the outer tube. In the heated section the couples are 
positioned In peripheral grooves in a hollow plastic insert 
which slips smoothly inside the core tube. The leads run 
down the center of this tube and out through a hollow copper 
tube which is soldered to the inside of the Inconel tube and 
serves as a current bus. In the unheated section the couples 
are similarly positioned but the leads are In axial grooves 
milled in the outer surface of the plastic insert. Details 
of the center tube construction are shown In Fig. III. A. 7- 
A similar correction for the temperature drop through the 
wall of the heated core tube is made. At the unheated end 
the plastic tube bearing the thermocouples Is removable to 
facilitate the installation of the center tube. 

The thermocouples were made from 30 gage, glass 
insulated, iron-constantan wire. Samples of the wire were 
calibrated against NBS thermometers. The maximum deviation 
from the standard tables 1 ^ was ± 0.25°F throughout the 
range of interest. 

The distribution of the thermocouple emf's to the 
potentiometer is through a multiple selector switch arrange- 
ment as shown in Fig. III. A. 8. The thermocouple leads at the 
test section are connected through MS type multi-pin connec- 
tors to copper extension leads from the switch console. The 
metal housing of the connector tends to Iso-thermalize the 
junctions to the copper. Details of the thermocouple circuit 

are shown In Fig. III. A. 9* 

The outside of the outer tube Is wrapped with about 

two inches of Fiberglas, then the entire assembly, Including 
the mounting channel is again wrapped with about one Inch of 
Fiberglas. The heat leak through the insulating blanket is 
determined by removing the core tube and allowing the system 
to come to equilibrium with no air flow. To minimize Intern- 
al free convection the tube is stuffed with Fiberglas. The 
analysis for the heat loss Is presented in appendix B con- 
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sidering the conduction loss down the current leads. This 
heat loss is actually only of significance when the outer 
tube is heated. 

The principal experimental interest is in the 
situation where only one of the boundary walls is heated. 

To minimize radiation interchange between the two surfaces 
a duplicate set of unheated tubes is provided in which the 
surface bounding the flow is nickel plated and polished. 

The derivation of the radiation interchange factors for the 
various geometries and the method of correction of the data 
for this loss are presented in appendix C. 

Additional details of the construction and operation 
of the apparatus are given by Leung ^ . 
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PIG. III. A. 2 - Power Supply and Instrumentation Circuit 
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PIG. III. A. 3 - Assembly Details 
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FIG. III. A. 4 - Air Metering System 



FIG. III. A. 5 - Outer Tube Thermocouple Locations 
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FIG. III. A. 7 - Core Tube Construction 








Switch Connectors 



PIG. III. A. 9 - Thermocouple Circuit 
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FIG. III. A. 10 Inconel Tubes Used in the 

Test Section 



FIG. III. A. 11 

Test Assembly, 
Exit End 
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Ill, b. Comparison of Experimental and Theoretical Values 


All of the geometries were run in laminar flow with the 
inner wall only heated. These results, together with the 
theoretical predictions, are shown in Fig. III.B.l. It is 
felt (see Appendix D) that these data have a probable un- 
certainty of + 5*1 per cent. All of the data are within this 
interval of the theoretical value and most of the points are 
considerably better. The tendency of the data to lie somewhat 
low at the first station (small x) may be due to an insufficient 
hydrodynamic entry length or to axial conduction in the tube 
wall. 

When the outer wall is heated the variation of the 
fundamental solutions with r* is much less pronounced. 
Accordingly, since the heat loss through the insulation in 
laminar flow was a significant fraction of the heat input 
only the smaller outer tube (l.O inches I.D.) was run in 
this condition. 

The results are shown in Fig. III.B.2 for r* = 0.376. 

The probable uncertainty in this data is + 7«3 per cent. 

Again all of the data lie well within this interval. It 
should be pointed out that the data for the outer wall 
heated do not really represent a constant heat flux situation 
but because of the dependence of the heat leakage on the 
tube wall temperature the heat flux declines about 10 per cent 
over the range of axial stations at which (®q Q ~ ^ mo ) was 
evaluated. The general good agreeement of the data with the 
theory suggests that If the axial heat flux variation is not 
too severe the local wall- to -mean temperature difference can 
be estimated from the constant flux solution without serious 
error. 

The apparatus Is so designed as to allow operation with 
simultaneous heating from both walls of the annulus. Data 
from this mode of operation are taken as examples of super- 
position of the fundamental solutions and are presented in 
appendix G. 



HEATED OUTER WALL HEATED 
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APPENDIX A 


CORRECTION FOR THE TUBE NALL TEMPERATURE 
DUE TO INTERNAL GENERATION 


As described in Section J jj. A, the thermocouples for 
measuring the tube wall temperature are In contact with the 
heated walls only on the side away from the flow boundary 
(i.e.j on the outside of the outer tube and the inside of 
the core tube). For the case where the wall is heated some 
difference between the inner and outer surface temperatures 
can occur. 

Consider a section of the tube as shown In Fig. A.l. 



FIG. A.l 

If we neglect temperature variation In the x-direction * , the 
equation of the temperature distribution in the tube wall Is 


The general 



solution 
t = 


dt , rq"' 

+ 3F + ~1T~ ~ 

to (a . l) is given by 

-n p 

r + Ginr + 


0 


H 


(A.l) 


(A. 2) 


where G and H are arbitrary constants 


*This same result is obtained if 


aft 

dx 2 


= 0, i.e 


lished flow with constant wall heat flux. 


fully estab- 
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For the heating of the outer tube we have the boundary 
conditions 


dt 

3r " 


0 at r 


r and t 
o 


t^ at r = r^. 


Using these we can determine G and H. Then, evaluating( A . 2) 
at r = r Q yields 


(t 


o 



At c 






!l!)](A.3) 


A somewhat more useful form of (A. 3) Tor the experimental work 
is 


At 

c 


1 

2?rk 




(A. 4) 


where (q/L) is the energy input per unit length of tube 
(BTU/hr.-ft . ) . 

When the core tube is heated the boundary conditions 

are 


dt 

'3F 


0 at r 


r. and t = t at r = r . 
i o o 


The application 
values of G and 


- t c'> * 


Alternatively 


At = oV 
c 2irk 


of these conditions yields the appropriate 
H in( A . 2) and we have 


n hi (r 

- ir [ — 


2 . r 2) 


- r i 2 ^] 


(A. 5) 



(A. 6) 
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For the geometries used in the experiments, the eval- 
uation of (A.4)and(A.6)yields the following correction formu- 
las : 

1 . Outer heated 

A. Inside diameter of outer tube 1.96 inches 

At = I .838 x 10 4 (q/L) (A*7) 

c 

B. Inside diameter of outer tube 1.00 inches 

At = 2.665 x 10 4 (q/L) (A. 8 ) 

c 

2. Core tube heated 

A. Outside diameter of core tube 0.50 inches 

At = 7.54 x lO -4 (q/L) (A. 9) 

c 

B. Outside diameter of core tube 0.376 inches 

At = 9.91 x 10 -4 (q/L) (A. 10) 

c 

In each case the wall temperature is given by 

4- = t (A. 11) 

u wall measured -At 

V-* 

It should be pointed out that for the laminar flow 
work presented here these corrections are entirely negligible, 
the range of (q/L) is between 10 and 50 BTU/hr.-ft. 
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APPENDIX B 


HEAT LEAK DETERMINATION FOR THE OUTER TUBE 

When the outer tube is heated the temperature difference 
between the wall and the surrounding ambient air causes 
some heat leak through the insulating blanket. For laminar 
flow, when the fluid film coefficient is relatively small, 
a significant fraction of the total energy generated passes 
out through the insulation or Is conducted down the power 
lead which is connected to the copper ring at the start of 
the heated section. This loss affects the experimental 
data In two ways, by altering the energy supplied to the 
fluid prior to any axial station at which we desire to com- 
pute the film coefficient, and by altering the local wall 
heat flux on which that film coefficient is based. 

To evaluate this loss, the Insulated outer tube was 
stuffed with FIberglas strips and, with the power on, al- 
lowed to come to equilibrium with the surrounding air. 

The power level was adjusted so that the temperature as- 
sumed by the tube was nearly equal to that which It would 
experience while running. The axial temperature distribu- 
tion was generally dome-shaped as shown in Fig. B.2, drop- 
ping off at the ends where the leads are connected and ris- 
ing to a maximum in the center. The dip in the center re- 
sults from the presence of a fourth Lucite holder which is 
not actually contacting the tube in this series of tests. 

For analysis, the system Is Idealized as shown In 
Fig. B.l. 


152 



Thermal Resistance 
of Air Film > 



FIG. B.l 


The temperature at any axial station is considered 
to be constant . (The results of Appendix A indicate that 
this is a good assumption.) The resistance of the insula- 
ting blanket and the free convection film resistance are 
lumped and considered to be a constant with length, that is, 
the heat flux going out through the insulation is assumed 
to be expressed by 


i t 


ins 


= h T (t - 


t ) 
a' 


(B.l) 


where the loss coefficient h^ is independent of temperature 
and axial position. The heat leak out the current lead and 
that into the unheated section of the tube are lumped and 
assumed to be 


^lead 


= Kr (t. 



(B.2) 
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where t^ Is the temperature at the connection of the lead to 
the tube . 

If we let t = t-t , the equation of the axial temper- 

cl 

ature distribution is 


a 2 t 

dx 2 


h L C t F 

kA 



0 


(B.3) 


where C^. = circumference of the tube, 

A = cross section area of the tube. 
The general solution to (B.3) is 


m A 


where 


t = 

2 

m 


G sinh(mx) + H cosh(mx) + r— p — 

Tt 

h L C t 
kA, 


(B.4) 


c 

If we take x = 0 when = 0, then G must be zero and we 
have 


t 


H cosh(mx) 


™ 1 A 


h L C t 


(B.5) 


By fitting (B-5) to the experimental curve at two points, 
h^ and H can be evaluated. We have the boundary condition 
that when x = x 

max 


" kA c 5x “ K L tm 

To obtain the lead leak constant, we differentiate (B.5) 
and substitute in (B.6), yielding 

mKA 

K t = -H sin h(mx ) * T . . - v \ 

L v max' (t) x = x • (B-7) 

max 

Typical values obtained for h T and K r are as follows: 
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Tube 
1.96 r ' 
1.00 " 


h L 

0.204 BTU/hr-ft 2 - °F 
0.284 


0.0493 BTU/hr-°F 
0.0381 


A comparison of the temperature distribution computed 
from (B.5) and the measured distribution is given in Fig. 
B.2. The value of t from (B. 5 ) Is forced to fit the exper- 
imental data at the start of the heated region and at the 
point of zero slope. 

When the core tube is heated the temperature of the 
outer wall is not greatly different from the ambient , but 
when the outer wall is heated t is about 50°F. For this 
condition the heat leak can be as much as 30# of the total 
energy input. 
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APPENDIX C 


THE EFFECT OF RADIATION INTERCHANGE 
BETWEEN THE TWO SURFACES 


If the inner and outer surfaces are each at a constant 
hut different temperature, the net heat exchange due to 
radiation is given by 


q 


r 


1 + 



1 ) 



(C.l) 


where the subscript i refers to the inner wall surface , 
o to the outer. 

During the tests, the axial temperature gradient 
is quite shallow, so that the constant temperature assump- 
tion Is not a serious handicap except very near the start 
of the heated length. The value of the emissivity of 
inconel was assumed to be 0.35*.> ° f the nickel-plated sur- 
face 0.1. Then, If we write (C.l) as 


q 


T1 

ri 



F 



(C.2) 


we can evaluate 


F 


e 



(C.3) 


for the various geometries. 

Note that the heat flux calculated from (C.2) is that 
at the inner wall, the flux at the outer wall is given by 


— 2(5 

*This is the value given by McAdams for a clean, 

well-polished surface, after repeated heating and cooling. 
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(C. 4 ) 


q 


ir 

ro 


r*q 


1 1 

ri 


r* 

Heated 

wall 


0.192 

outer 

O.O965 


inner 

0.218 

0.256 

outer 

Inner 

0.0953 

0.1932 

0.376 

outer 

inner 

0.0933 

0.160 

0.501 

outer 

Inner 

0.0915 

0.1352 


The objective of the adiabatic wall tests is to de- 
termine one of the fundamental solutions directly. The 
effect of this radiation loss on this determination is 
shown by the superposition example of I.B. 

1 • Inner wall heated 

The experiment measures t, and determines t , but 

l m 

from (l.B. 17 ) and (I.B.18) we get that* 



V 


q o ( 0 lo ~ 0 mo) ^(C.5) 


It is ( 0 11 - ) that we seek. Here q^° = r*q ri " - q 


ii 

and q . M = q fl 


mi 

measured -0 ^! 


^ii - % 


x <*i - *m> 

} ^ "Vk 


r* q 


- 0 ) 
" V 1 0 m r\ * 


mo- 


(C.6) 


In laminar flow the temperature profile develops rather 
slowly, and for most of the range of axial distance covered 


* The experiment deals only with fundamental solutions of 
the second kind, the superscript (2) has been omitted for 
brevity. 
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In the experiment , e ±Q is nearly zero. For the purposes 
of the correction to the experimental data it Is assumed to 
he so. Note that this is making the correction to the mean 
temperature due to heating at the outer wall, but assuming 
that the heat supplied at the outer wall does not affect 
the Inner wall temperature 
2. Outer wall heated 

In an analogous manner, (I.B.17) and(l.B.19) can be 
combined to give 


^0 ■ 


> I Vk 




where = q" measured + r *<j" rl - <i£- 

Equation (C.7) Is used to correct the outer wall heated data 
In the same manner as (C.6), i.e., is assumed to be zero 

For laminar flow. Inner wall heated (the worst case) 
the heat transfer due to radiation can be as great as 10$ 
of the total heat transferred. 
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APPENDIX D 


EXPERIMENTAL UNCERTAINTY 


The experimental uncertainty in the fundamental solu- 
tions, 9 - 0 , and the position parameter x, can he esti- 

ill ? -y 

mated using the method described by Kline and McClintock' . 

1 . Uncertainty in x~ 

The defining equation is* 


x 


v\ 

RePr 


(D.l) 


The uncertainty in is negligible. The axial position of 

the thermocouples is easily accurate to ± 1/16 inch. This is 
only ± 2 % of the first station evaluated, and becomes negli- 
gible as x grows larger. The uncertainty in Prandtl number 
Is about ± 2.55*. The Reynolds number is evaluated from 

Re - (]£> TT- < D - 2 > 

The determination of the flow rate, w, is subject to the man- 
ometer error and the uncertainty of the orifice coefficient. 

The orifice coefficient should be reliable to + 0.55* (see 
Ref. 23). Even though the pressure differentials are small 
(0.5 inch) the manometer uncertainty is of the order of + 1 % 
or ± 0.55* when entered as the square root. The uncertainty of 
the density of air Is about + 1 %. These combine to create a 
probable uncertainty in the flow rate of + 1.25*. The viscos- 
ity of air is known to + 1.55** hence the probable uncertainty 
in Reynolds number is + 1.75*. The uncertainty in 5c becomes 
+ 3 . 65 * at the first station, declining to + 3 . 05 * down the tube. 


*Note that x = 0 corresponds to the start of the heated section. 



2. Uncertainty in 9 ~ e m 

The computing equations for 9 - are (C.6) and(C.7)* 

a . Heating at the inner wall 

During the tests, t ± - t m was set at approximately 

30 °F, The uncertainty in t^ is probably -0.5° • t m is computed 

by adding the temperature rise due to heating to the inlet 

temperature. This temperature rise reaches a maximum of 20°. 

The ammeter and the voltmeter used for the power measurements 

have operating uncertainties of iO.75# and ^2# respectively. 

The variation in wall thickness is ±0.75#* The uncertainty 

in power measurement is then ±2.3#.> and the uncertainty in 

t about tl°. The relative uncertainty in t. - t becomes 
m ■*" ^ 

about ±3*5#. The radiation loss is about 6# of the total 

input; if this estimate is off by -25#, the error in q_" is 

+1.5#. The correction term in (C.6) is, at its maximum, only 

about 3# of the first term, so a ±50# error here means only 

±1.5#. The thermal conductivity of air is probably reliable 

to ±2#. The combination of these effects produces a total 

probable uncertainty in 9 ^ of ±5 -l#* 

b . Heating at the outer wall 

The same considerations apply to 0 QO - except 

that the heat leak through the insulating blanket is 20# to 
30# of the total input. It is felt that the estimate of 
heat leak is good to ±10#, hence the uncertainty in wall heat 
flux becomes ±3-8#. This makes the uncertainty in t Q - t m 
±5.6#, and the probable uncertainty in 9 qq - 0 mo ±7*3#* 

A summary of the uncertainties is given below. These 
uncertainties in all cases are based on 20 to 1 odds. 
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Item 


Core tube heated 


Thermal conductivity 

Viscosity 

Prandtl number 

Flow rate 

Re 

x 

Heat flux 


- 2 . 0 $ 

* 1 . 5 $ 

* 2 . 5 $ 

* 1 . 2 $ 

*1.9$ 

*3$ to *3-6$ 
* 2 . 8 $ 
^ 5 - 1 $ 


Outer tube heated 


* 2 . 0 $ 

*1.5$ 

*2.5$ 

* 1 . 2 $ 

*1.9$ 

*3.0 to *3.6$ 
*3.8$ 

*7-3$ 
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APPENDIX E 


DETAILS OF THE NUMERICAL SOLUTION 


The method of Berry andde Prlma^ Is a convergent itera- 
tion scheme for finding the proper eigenvalues for a Sturm- 
Llouville equation. For the case at hand, this equation is 
given by (II.F.2). Actually, (II.F.2) represents a system 
of equations since two different boundary conditions are ap- 
plied for each of the fundamental cases, and all the condi- 
tions are repeated for each value of radius ratio considered. 

For each step of the Iteration we seek a solution to 
( II.F.2) subject to (iI.H.l). The requirement that this so- 
lution satisfy the boundary condition only at r = r*, creates 
an initial value problem. With ^ = ( A n ) k a value is 

chosen for R n (r*) or R n (r*) (note that only one of these is 
non-zero for any given case) and a numerical Integration of 
II.F.2 is performed from r*<r<l . In general, the solution so 
generated will not be normalized, that Is 

I T l w(F) [R n (r) ] 2 dr j 1 
Since Eq. ( II.F.2) is linear, it follows that 


^ R n^))k= Vn^) 


where 


K 


n 


/ r *w(r) [R n (r) ] 2 d p 


(E.l) 


(E.2) 


will be the desired solution satisfying (II.H.l). Equation 

2 

(II.H.2) supplies the (k + l)st approximation to A n and the 
iteration procedure continues. The values given by the 
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asymptotic expressions (see Section II. G) are used as the 

2 

first approximations to A^ in all cases. The correction 
supplied by (II.H.2) is very strong; after the third iteration, 
the residual, given by (R ( l) ), (R' (l) ), , is less than 

il K 11 .TV ^ 

10.0005 which is used as the convergence criterion for A . 

The actual numerical integration is performed using a 

24 

three-point predictor-corrector method due to Hamming . To 
illustrate the method, refer to Fig. E.l. 


01 2 3 (m-l) m 

r = r* r = 1 

FIG. E.l 

The interval r*< r < 1 is divided into equal increments 
Ar. The station r = r* is numbered 0, r = r* + Ar is numbered 
1, etc., up to r = 1, which is station m. Omitting the iter- 
ative index, subscript k, let the value of nth eigenfunction 
at station p, p = 0, l,2,...m, be denoted R n (p) and the deriv- 
atives by the prime notation R^p) and R*(p) . After the 
scheme is started, we predict the derivative at the p + 1 
station by 

R^p + 1) - R^P - 3) + |(AF)[2R^' (p) + 2R^(p - 2) 

-R"(p - 1) ] (E.3) 

Using this we can calculate 

R n (P + 1) = jj/9R n (p) - R n (p - 2) + 3(Ar) [R^(p + l) 

+2R^(p) - R^(P - l)]j (E.4) 
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The value of R^(p + l) can be calculated from the differential 
equation and the corrector (E.4) applied to both R^(p + l) 
and R n (p + l)- These next better values are then used to 

recalculate R"(p + l) and(E.4) applied again to settle the 
values of R^(p + l) and R n (p + l). The scheme goes as fol- 
lows : 

(1) Compute R^(p + l) from (E.3) 

(2) Compute R^(p + l) from (E.4) 

(3) Compute R^(p + l) from (II.F.2) 

(4) Compute R^(p + l) from (E.4) (adding another 

prime to every term) 

(5) Compute R n (p + l) from (E.4) 

(6) Repeat steps 3* 4, and 5* 

PR 

Equation (E.3) is the predictor used with the Milne 3 method. 

Considerable care is used in starting the procedure and 
an iteration is performed to correct the first station values 
by 

R n (l) = V 0) + I5R^(0) + 8 r^(i) - R^(2)] (E.5) 

After the integration is complete, the norm of the 
eigenfunction 

(i w(F)[R n (F)] 2 dF 

is calculated using Simpson's rule. When the iteration 
scheme has satisfied the criterion for closure, the eigen- 
constants are calculated from (II.F.7)* again using Simp- 
son's rule for the integration. 

To provide a validity check on the method, the first 
six eigenvalues for the circular tube with constant wall 
heat flux (the fundamental solutions of the second kind for 
r* = 0) were calculated. Table E.l is a comparison of 

these values with those published by Siegel, Sparrow, and 
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20 

Hallman . These integrations were performed using 250 
divisions of the interval. 


n 

Siegel, et al 

Present Values 

1 

7.16654 

7.16654 

2 

12.9508 

12.9508 

3 

18.6638 

18.6636 

4 

24.3530 

24.3532 

5 

30.0317 

30 . 0317 

6 

35-7040 

i 

35.7030 


The machine program is written in the Burroughs 1 ver- 
sion of ALGOL. The following symbols are used in the pro- 
gram. 


Quantity 

Machine Symbol 

A 

L 

n 

A 2 

LL 

n 

Weight function, w( 

» W(I) 

Position index, p 

I 

R n (p) 

R(I) 

R’(p) 

RP(I) 

R>) 

RDP(l) 


Ar 

DEL 

sf* 

ro 

C 
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COMMENT LAMINAR EIGENVALUES FOR THE ANNULUS; 

INTEGER N,I; 

ARRAY RBR(3^0),W(340),G(3^0),R(340), RBRSQ(340), RBRHF(340), 
LNR(34o), FRl(340) , FRO(3^0), RP( 340 ) ,RDP( 340 ) ; 
INPUT VALUES ( L , RO , RPO , DEL , EPS 1 , EPS 2 , THETA I , THETA 0 , RSTAR , K ) ; 
OUTPUT ANSWR (L,LL,RO,RPO,CI,CO,RSTAR,FOR I=(l,l,N) ; ( (C*0 . 5) • 

R(i) ) ) > 

NOTE (L), CHEK1 (R(n) ,PR(N) ) , CEK2 (NORM); 

FORMAT FRMT(B5,*L*,X12.8,B2,*LL*,X12.7,B2,*R0*,X12.8,B2,*RP0*, 
X12.8,B2,*CI*,F15.8,B2,*C0*,F15.8,B2,*RSTAR*,X6.3,W2,125 
(8X14 . 10,W2 ) ) y FRMT2 (*THIS INTERMEDIATE L IS* ,Xl4 . 7 ,W2) ; 
FORMAT FRMT3(*R(N)*,X15.10,*RNP*,X15.10,WO) , 
FRMT4(*NORM*,X15.10,WO); 

PROCEDURE CONST(M,K,DEL,W( ) ,F( ) ,R( ) ;C) ; 

BEGIN INTEGER I,M; C=Q.O; FOR I=(l,2,M-2); C=C+DEL(4W(l) .F(l) . 
R(I)+2W(I+1).F(I+1 ) .r(i+i))/ 3; C=(C+4DEL.W(M-l).F(M-l) . 
R(M-l)) (k*0.5)/3j 

RETURN END; 

COMMENCE . . READ ( ; ; VALUES ) ; 

DR=1 . 0 ; 

IF (DR NEQ RSTAR); BEGIN DR=RSTAR; 

LGRS=LOG( RSTAR) ; B=( (RSTAR*2)-l)/LGRS; 

EN=l/( 2 ( l+( RSTAR* 2 ) -B ) ( ( 1 -RSTAR ) *2 ) ) ; 

M=(l -RSTAR) /DEL; N=FIX(M); 

FOR I=(l,l,N-l); BEGIN 
RBR ( I ) =RSTAR+I . DEL ; 

RBRSW ( I ) =RBR ( I ) .RBR(l); 

RBRHF( I ) =SQRT(RBR( I ) ) ; 

LNR ( I ) =LOG ( RBR+I ) ) ; 

w(i)=en(rbr(i)-(rbr+i)*3)+b.rbr(i)lnr(i)) ; 

fri(i)=(-4rstar.en/(i+rstar))(( (i-b)/4)(rbr(i)*2-lnr(i))-(rbr 
( I ) * 4 ) /l6+B( RBR ( I ) *2LNR ( I ) /4- ( RSTAR*2 ) ( l+( RSTAR*2 ) -B) 

/4+ ( 1+LGRS ) /4 - ( RSTAR*4 ) /l6 ) -THETA I ; 
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fro(i)=(-4en/(i+rstar)J ( ( ( l-B)A) (RER(I)*2-1 )-(RBR(i)*4-1)/i6 

+B(RBR(l)*2)LNR(l) /4-(RSTAR*2) (RSTAR*2-B)LNR( I )/4 ) -THETA 0 
END END.;W(N)=0.0; RBR(N)=1.0; RBRSQ(N)=1 .0; RBRHF(N)=1 .0 ; 

INF. .LL=L. L; 

FOR I=(l,l,N); G( l)=LL.W( i) ; 

RDP0=0 . 0 ; RP( 1 )=0 . 0; R(l)=ROj 
FOR P=l. 0,2.0; BEGIN 

RDP( 1 ) =- (rp( i )+g( 1 ) .r( 1 ) )/rbr( 1 ) ; 

RP(l)=DEL(RDPO+RDP(l) )/2; 

R ( 1 ) =RO+DEL ( RPO+RP ( 1 ) ) /2- ( DEL* 2 ) ( RDP ( 1 ) -RDPO ) /l 2 END ; 

RP( 2)=R/ ( 1 ) =DEL( 3RDP( 1 ) -RDPO)/ 2; 

R(2)=R0+(RP0+4RP(l)+RP(2) ) DEL/3; 

FOR P=l. 0,2.0; BEGIN 

RDP(2)=-(RP(2)+G(2) .R(2))/RBR(2)j 

RP ( 1 ) = ( 5RDP0+8RDP ( 1 ) -RDP ( 2 ) ) DEL/l 2+RPO ; 

R( 1 ) = ( 5RPO+8RP ( 1 ) -RP ( 2 ) ) DEL/l 2+RO; 
RDP(i)=-(RP(i)+G(1).R(1))/RER(1); 

RP ( 2 ) =RPO+ ( RDP0+4RDP ( 1 ) +RDP ( 2j ) DEL/3 ; 

R( 2)=R( 1 )+DEL(RP( 2)+RP( 1 ) )/ 2- (DEL* 2) (RDP( 2)-RDP( 1 ) )/l2 END; 

RP ( 3 ) =RP ( 2 ) +DEL ( 3RDP ( 2 ) -RDP ( 2 ) -RDP ( 1 ) ) /2 ; 

R(3)=(9R(2)-R0+3DEL. (RP(3)+2RP(2)-RP(l)))/8; 

FOR P=l. 0,2.0; BEGIN 

RDP(3)=1(RP(3)+G(3).R(3))/RBR(3); 

RP(3)=(9RP(2)-RP0+3DEL.(RDP(3)+2RDP(2)-RDP(l)))/8; 

R(3)=(9R(2)-R0+3DEL.(RP(3)+2RP(2)-RP(i)))/8 END; 

RP(4)=RP0+(4.0)DEL.((2.0)(RDP(3)+RDP(l))-RDP(2))/3.0; 

R(4)=(9R(3)-R(l)+3DEL.(RP(4)+2RP(3)-RP(2)))/8; 

FOR P-1. 0,2.0; BEGIN 
RDP(4)=-(RP(4)+G(4).R(4))/RBR(4); 

RP(4)=(9RP(3)-RP(1)+3DEL. (RDP(4)+2RDP(3)-RDP(2)))/8; 
R(4)=(9R(3)-R(i)+3DEL.(RP(4)+2RP(3)-RP(2)))/8 END; 

FOR I=(4,1,N-1); BEGIN 

RP( I+1)=RP( 1-3 ) + ( 1 . 3333331 ) DEL. ( ( 2.0) (RDP(I)+RDP( 1-2) )-RDP( 1-1 ) ) ; 
R(l+l)=(9R(l)-R(l-2)+3DEL. (RP(l+l)+2RP(l)-RP(l-l) ) )/8; 
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FOR P=l. 0,2.0; BEGIN 

RDP(I+1)=-(RP(I+1)+G(I+1) .R( I+l) )/RBR(l+l) j 

RP(l+l)=(9RP(l)-HP(l-2)+3DEL.(RDP(l+l)+2RDP(l)-RDP(l-l)))/8; 

R(l+l)=(9R(l)-R(l-2)+3DEL.RP(l+l)+2RP(l)-RP(l-l)))/8 end END; 
WRITE ( ; ; CHEK1 , FRMT3 ) ; 

N0RM=0.0; 

FOR I=(l,2,N-l); NORM=NORM+DEL( 4W( i) (R( I) -R( i) )+2W( I+l) (R( I+l 
).R(l+l)))/3; 

WRITE ( ;;CEK2,FRMT4); 

c=i/norm 

IF ((ABS(C.R(N).RP(N))) GEQ, EPS3)j GO RECORD; 

CONST (N,C ,DEL,W( ) ,FRl( ) ,R( ) jCl) ; 

CONST (N jC ,DEL,W( ) ,FRO( ) , ,R( ) ;CO) ; 

WRITE ( ; ; ANSWR , FRMT ) ; GO COMMENCE; 

RECORD. .WRITE ( ; ;N0TE,FRMT2) ; 

L=SQRT(LL+C.R(N) .RP(N)); 

GO INF; 

FINISH; 
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APPENDIX F 


DERIVATION OF THE EXPANSION COEFFICIENTS 

In order to establish completely the solution to 
the energy equation it was necessary to expand the initial 
condition on 0^ ^ as an infinite series of the eigenfunctions 
(II. F. 6) with the expansion coefficients, C , given by 
(II.F.7). For computational purposes, particularly to 
establish the asymptotic values of the C^ , it was necessary 
to provide a more convenient form for the integrals appearing 
in (II. F. 7). 

We desire to evaluate the Integral 

1 P _ 

/ w(r) R^(r)dr (F.l) 

r>* n 


which Is the norm of the eigenfunction. Let us write 
Eq.(lI.F,2) as 

d _ dR o „ 

— (r ) + h w(r) R = 0 (F.2) 

dr dr n 

Taking of the above we have 


3 

ax 


a aR 

— (r^) 

L ar Sr J 


+ 2A w(r) R n + w(r) 


aR 

n 

ax 


0 


Since the order of partial differentiation may be reversed 
this may be written as 


a - 
a? _ 



aR ■ 

(r- 2 ), 

ax J 


_ p _ aR 

+ 2X w(r)R + w(r) — - 

n ax 


0 


(f . 3 ) 
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If we multiply ( F . 3) by R n and integrate between r* 
we have 


id d dR _ l_p_ 

/ R — tr — ( — -)] dr + 2A / w(r) R„ dr 
r* n dr dr dA r* 


_ 1 dR 

+ A 2 / w(r) ( — -)R dr = 0 
r* dA 


Integrating by parts twice we can get 


1 o d dR x 

2A / w(r) T dr = -[ r R (r) ( — “) ^ 

r* n n dr dA r * 

dR dR 1 1 dR r d dR 

+ [ (-*) r (-^) 1 - / (-^) [ — (v -?) 


dA 
2 /-> 


dr 


>* dA dr dr 


+ A^ w(r) R ] dr = 0 


If A = A the Integral on the right hand side of 
vanishes and we have* 


rrr\ t -.2 


n 

2A 


n 

dR 

dR 

n } 

T (— 

dA 

dr 

n 



d dR 
- t r R (r) — (— ) ] 
n dr dA 


n r’ 


there the derivative dR^/dA^ is used to indicate 

(aV^-x • 


n 


and 1, 
(P.^) 


(p .5) 
.5) 

(F.6) 
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For the various fundamental cases certain of the terms in 
(F.6) will vanish by virtue of the boundary conditions but 
for each specific case the value of the norm may be readily 
obtained. 

The numerator of (II.F.7) has the form 


1 


/ w(r) 
r* 


(r) R (r) dr 
j ' n ' 


(F.7) 


where Fj k ^(r) = 0j k ^(r,o). To evaluate this, multiply 
(F.2) by F(r) and integrate by parts. Then 


J 1 P^(r) 

-p* J 


w(r) R n (r) dr = 



1 SR dF^ k ) 

- / (r-^) (-4—) dr 

r* Sr dr 


/, \ ^R 

p k) (?)(4) 

J Sr 


1 

] 

r* 


(F.8) 


Let us consider (F.8) for each of the initial conditions 
developed in section II. D. 

For the first case 


^i 1 ) (r,o) = F^tF) 


in r 
in r* 


, then 


dFp-^r) 

dr 


1 

— so that 

r in r* 





) w(r) R n (r) dr 



(- 4 ) 

Sr — 
r= 


* 


(F . 9) 
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Given that F^^r) = - 1 we obtain in an analogous 

manner 


f F^Ur) w(r) R (r) 

jt O A i 

r* 



dR 

(-d 1 ) 

7=1 


(F.10) 


For case 2, the first term on the right hand side of (F.8) 
vanishes by virtue of the boundary conditions and we have 




( 7 ) 


w(r) R n (r)dr 


1 

I r 


n 


dR dF^ 2 ) 

(_J1) (_J_) dr 

dr dr 


Integrating by parts yields 


/ f( 2 ^ (r) w(r) R (r) dr 



dF 


( 2 ) 


7 (-4—) V r ) 1 

dr 


1 1 d dF^ 

— o / R n (r) — (r —± ) dr 

A r* dr dr 

n 


Recall that F^(r) = 


- d[ 2 ^ (r ) )** and 


J 


fd 


,(2) 


,( 2 ) 


4 (7 “Tv, . ( ^i_) „ ( f) 

dr dr dx 


so that we have 



1 ,(2) 1 dF 


( 2 ) 


/ f!j (r) w(r) R (r) dr = — p [ r (— J — )R (r) 1 

-¥• J II II 


n 


dr 


i i de 


( 2 ) 


+ 


/ ( — r 11 -) w(r) R (r) dr 


A r* dx 
n 


n 


Further 


w(r) R n (r) 


1 d dR 

— (r -~ ) 

dr 


A 2 dr 


n 


so that 


1 (2) - - _ 1 dF^ 2) 

/ Fj (r) w(r) R n (r) dr = + — g [r(— — ) R n (r)] 


V dr 


i de< 2 ) 


dR 


A dx 

n 


'w' ^ J W 

(_aJ_) [ r _n ] 

dr 


The second term on the right vanishes and 


dF*; 2 ) 1 

— jL — = + 

dr 2(l-r*) 


at r = Ty and zero on the 
opposite wall. 


Hence 




(r) w(r) R n ( r) dr = 


r* R (r*) 
n ' 

A 2 2(l-r*) 
n v ' 


(F.ll) 
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and 


/ F^ 2 )(r) w(r) R n (r) dr = - 


R (l) 
n 7 


A 2 2(l-r*) 
n 


dF 


(3) 


! o ) _ ^ -J 

For the third kind F^'(r) = -1, then — ± — = 


and we have immediately 


dr 


/ F$ 3 )(r) w(r ) R (r) dr = - 


r* SR 

— 2 (-^L 

1 dr r=r* 


n 


and 


1 SR. 


f 1 F^ 3) (r) w(r) R (r) dr = — ( 7 =?)- 

jt U Al ^ 


A dr r=l 
n 


For the fourth kind 


F (4) (7)=-i £nr 

^ 2(l-r*) 


and 


F 


n^)( r ) = - 

o ' 


r 

in — 


2 ( 1-r*) r^ 


so that (see case l) 


(F.12) 

0 


(P.13) 


(F.l4) 
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(P.15) 


r 


f 1 F$ 4 )(7) w(r) R (r) dr = 
* 1 “ 


r* R r ( r*) 
\ 2 2( 1-r*) 


and 


/ p( 4 )(r) w(r) R n (r) dr = 

r* 0 


R (l) 
n 

\ 2 2(l-r*) 
n 


(f.i6) 
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APPENDIX G 

EXAMPLES OP SUPERPOSITION OP THE FUNDAMENTAL SOLUTIONS 


The graphs of the fundamental solutions. Pigs. II. I. 1 
to II. I. 12, are particularly convenient for rapid calculations 
involving superposition of the elementary boundary conditions. 
To illustrate this computation let us consider some specific 
examples. 


1. Heat flux specified on one wall, temperature on the other 

Assume that the heat flux on the Inner wall Is as shown 
in Pig. G«l, while the outer wall is held at 212 °P. The fluid 
is air which enters at 70 °P with a Reynolds Number of 1428. 

The outer tube is 4 inches In Inside diameter, the core tube 
1 inch in outside diameter. The tubes are 150 inches long. 

We seek the temperature at the Inner wall, the fluid mean 
temperature, and the heat flux at the outer wall as functions 
of axial position. The pertinent physical constants are: 

r* = 0.25 

D ■ 0.25 ft 
h 

RePr = 1000. 

This problem will require the superposition of cases 
3 and 4. Let us approximate the cosine distribution by the 
steps shown In G.l. If £ is the value of x at which 
the step of heat flux, q^, Is applied, we will have the 
following; 


n 

W 

0.005 

0.015 

0.035 

0.045 

0.050 


<n V* 
10 
100 
90 

- 90 
-100 

- 10 


- °F 
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We can obtain the inner wall temperature from 


t t ra = 70 + 142 0<3>(x) + 2( q” n D h /k) /^(x - e n ) 


(G.l) 

The summation is taken only of those steps for which x > £ . 
Similarly, the outer wall heat flux is given by 


k 


q"(x) = T 3 - t 142 ®^<3+ S(q^ ni ji 
h 


) ®^ ) (x 


4 n )] (G.2) 


and the mean temperature by 


t m (x) = 70 + 142 e^(x) + - £ r ) (0.3) 

These values are shown in Pigs. G.2 and G.3. Examples 
of the calculation are as follows. 


1 .) 


x = 0.01 

t ± = 70 + (l42)(0) + (10) (0.106) + (100) (0.0868) 
t ± = 79.7 °P 

q^ = (o. 0146/0. 25)[ (142) (4.28) - (10) (0.0004)-(l00) (0) ] 

q" =609 Btu/hr ft 2 
o 

t m = 70 + (142) (0.225) + (10) (0.008) (100) (0.004) 
t = 102. 4 °f 
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2.) x = 0.05 

t ± = 70 + ( 142) (0.313) + (10) (0.165) + (100) (0.l6l) 

+ (90)(0.l49) - ( 90) (0.119 ) - (lOO) (0.87 ) 
t ± = 126.3 °F 

q£ = (0. 0146/0. 25)[(l42)(l. 85) - (l0)(0.0782) - (100) 

(0.0670) - (90)(0.0420) + (90)(0.0034) + (lOO) (O)] 
= 251 Btu/hr ft 2 

t = 70 + (i42) ( 0.568) + ( 10) (0.035) + ( 100) (0.0328) 

+ (90) (0.0265) - (90) (0.012) - (100)(0.004) 

t m = 155.2 °F. 

Notice that in this example, even though heating Is 
taking place at the Inner wall, the Inner wall temperature 
Is below the mean temperature. This Is due to the strong 
influence of the hot outer wall. 
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BTU/hr - ft 








2. Heat flux specified on both walls 

This is a situation which could he set up experimentally. 
The inner wall heat flux is essentially constant at an 
arbitrary level but the flux at the outer wall decreases 
with x since the heat loss through the insulation increases 
with the wall temperature. 

The radius ratio is 0.376. The values of the funda- 
mental solutions are given in Pigs. G.8 and G.9. These 
values were taken from cross plots made up from Tables II. 1.5 
and II. 1.6. 

The results of two asymmetrically heated runs are shown 
in Pigs. G.4 through G.7. In run G.l the wall heat fluxes 
were approximately equal, in run G.2 the heat flux at the 
inner wall was approximately twice that at the outer wall. 

The actual wall fluxes as they were measured (corrected for 
radiation interchange and heat leak) are shown in Pigs. G.4 
and G.6 together with the approximate distributions chosen 
to represent them. A comparison of the actual and predicted 
temperatures is shown in Pigs. G.5 and G.7. 

Por run G.l, the fluid mean temperature is computed 

from 

t m (x) = 75.7 + 178 eJjXx) + g(q" V k )n {0A) 

where £ n is the value of x for which the n-th step of 
heat flux is applied. Similarly 

t x (x) - t m (x) + I78 te^(x) - e^XxJl 

+ £ K <*.«»)] (0.5) 
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and 


t o ra - t m (x) + g(qj V k )n [ 9 io )(x -y - 


+ 178[0^^(x) - e^)(x)] 


( 0 . 6 ) 


The summations Include only those steps for which x > t n « 
For run G.2 the mean temperature Is given by 


t m (x) =82+2 (q£ V k >n (x ' e m> 


+ g(lS V k > 0 io’< x " e no) 


(G.7) 


The wall temperatures are given by 


tl (x) = t m (x) + g(qj V k ) n [8 W < X - 5 ni> - ®ii )(x - W J 


S K V k >n [e i^ (x - 4 no> 


«no> ] 


(G.8) 


t 0 Cx) - t m (x) + 2 V^)n [e oi )( * - 4 no> ' B lo ' «no> 1 


+ 2 (q- V k )n t8 oi ) ( x - *ni> ’ 6 if < x 


^nl^ ^ 


(0.9) 
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As an example of the calculation take the location 
x = O. 07 . 

t m = 82 + 200 ( 0 . 075 ) - 13(0.0109) + 105(0.204) - 15(0.102) 

t = 116.7 °F 
m 

t ± = 116.7 + 200(0.147) - 13(0.112) - 105(0.091) + 15(0.0759) 

t ± = 136. 2°F 

t Q = 116.7 + 105(0.202) - 15(0.190) - 200(0.0335 + 13(0.011) 
t = 128.4 °F 

Notice that the terms for which > x are all zero. 

As a comparison the wall temperatures in run G.l 
were computed assuming the wall fluxes were both constant 
at their average values. This prediction Is shown as the 
dashed line In Fig. G.5. While this prediction represents 
the wall temperatures less well than the multi-step approxi- 
mation. It is still reasonable since the actual axial variation 
of the heat flux is not severe. Notice that the constant 
flux prediction deviates from the measured values most on 
the outer wall, near the step, where the actual wall flux 
is varying most rapidly. 

The wall temperatures in run G .2 were computed assuming 
the heat fluxes to be those which were actually measured but 
that the wall-to-mean temperature differences were those 
given by the fully developed (large x) solutions. This 
method provides an adequate prediction for the wall tempera- 
tures for x > 0.1 but, again, the axial variation of flux 
is not very severe. 
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For the fundamental solutions of the second kind, the 
first significant eigenvalue is larger than for the other 
cases = 0) , and, hence, fully developed conditions 

are achieved most rapidly for this case. 
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Constant Flux Prediction 



FIG. G.5 
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The Fundamental Solutions of the Second Kind 

r* = O.375 



FIG. G.8 
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The Fundamental Solutions of the Second Kind 

r* = 0.375 
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